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A  variational  development  of  the  governing  equations  and  boundary  conditions  for  the 
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for  each  of  the  sandwich  facing  elements,  in  terms  of  either  a  generalized  force  or 
displacement,  a  condition  for  either  the  shear  loading  or  the  average  displacement  of 
the  core  Is  obtained.  The  governing  equations,  specialized  for  small-displacement 
buckling,  are  applied  to  the  case  of  axial  compression  of  circular  cylindrical  shells 
The  unified  theory  analysis  reflects  three  modes  of  buckling:  classical  general  in¬ 
stability  and  both  symnctric  and  antisymmetric  wrinkling.  Comparison  of  the  results 
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in  good  approximation  by  a  linear  theory. 
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SUMMARY 


A  variational  development  of  the  governing  equations  and  boundary  condi¬ 
tions  for  the  bending  and  buckling  of  sandwich  curved  plates  and  shells 
is  presented.  The  theory  is  applicable  to  sandwich  construction  of  the 
honeycomb  type  and  may  be  used  to  describe  bending,  buckling,  and  post- 
buckling  behavior.  With  the  coupled  effects  of  core  transverse  shear 
and  normal  strain  included,  the  theory  reflects  the  presence  of  8  degrees 
of  freedom  and  2  constraining  conditions.  Nine  boundary  conditions  per 
curved  plate  or  shell  edge  can  be  specified;  apart  from  the  four  of  the 
usual  type  for  each  of  the  sandwich  facing  elements,  in  terms  of  either 
a  generalized  force  or  displacement,  a  condition  for  either  the  shear 
loading  or  the  average  displacement  of  the  core  is  obtained.  The  gov¬ 
erning  equations,  specialized  for  small-displacement  buckling,  are  applied 
to  the  case  of  axial  compression  of  circular  cylindrical  shells.  The 
unified  theory  analysis  reflects  three  modes  of  buckling:  classical 
general  instability  and  both  symmetric  and  antisymmetric  wrinkling. 
Comparison  of  the  results  with  published  experimental  data  for  honeycomb- 
core  sandwich  shells  of  practical  dimensions  yields  very  strong  evidence 
to  suggest  that  the  axial  compression  buckling  load  and  failure  mode  of 
such  shells,  unlike  homogeneous,  isotropic  shells,  are  given  in  good  ap¬ 
proximation  by  a  linear  theory. 


FOREWORD 


The  work  reported  herein  constitutes  a  portion  of  a  continuing  ef¬ 
fort  being  undertaken  at  Stanford  University  for  the  U.S.  Army 

$ 

Aviation  Materiel  Laboratories  under  Contract  DA  44-177“AMC-115 (T) 
(Task  1F162204A17002)  to  establish  accurate  theoretical  prediction 
capability  for  the  static  and  dynamic  behavior  of  aircraft  structural 
components  using  both  conventional  and  unconventional  materials.  A 
predecessor  contract  supported  investigations  which  led,  in  part,  to 
the  results  presented  in  References  2  and  4. 


^Redesignated  Eustis  Directorate,  U.S.  Army  Air  Mobility  Research 
and  Development  Laboratory. 
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INTRODUCTION 


Although  substantial  interest  in  sandwich  construction  for  aeronautical 
applications  was  evidenced  during  and  just  subsequent  to  World  War  II, 
it  has  been  only  in  the  last  few  years  that  the  inherent  potential  of 
composite  construction,  in  general,  has  been  truly  appreciated.  The 
continual  stri-ving  for  lightweight  structures  with  higher  strength  and 
stiffness  has  caused  significant  interest  in  composites  for  aerospace 
structures  applications  and,  as  a  result,  has  led  to  renewed  interest 
in  the  sandwich  type  of  composite.  More  and  more  aerospace  designs  are 
reflecting  the  usage  of  sandwich  construction  in  primary  structural 
members  as  fabrication  and  prediction  capabilities  improve. 

As  is  usually  the  case  in  aerospace  developments,  engineering  achieve¬ 
ments  outstrip  theoretical  prediction  of  performance.  Thus,  to  avoid 
costly  experimental  design  and  conservative  analysis  techniques  and  to 
realize  the  potential  of  composite  structures  in  competition  with  con¬ 
ventional  metal  structures,  theory  and  practice  must  be  made  concurrent. 
Unfortunately,  the  price  to  be  paid  for  achieving  more  optimum  structures 
is  complexity.  Since  the  inherent  complexity  manifests  itself  in  the 
fabrication,  testing,  and  analysis  of  sandwich  construction,  each  of 
the  areas  must  receive  comparable  attention  and  the  problems  must  be 
attacked  directly  rather  than  accommodated  by  the  introduction  of  con¬ 
servative  practices.  The  status  of  the  design  criteria  for  stiffened 
metal  shells  attests  to  the  penalties  incurred  by  the  absence  of  a 
satisfactory  theoretical  solution  of  the  problem,  even  though  the 
fabrication  and  testing  of  such  shells  is  very  well  in  hand. 

A  perusal  of  Reference  1,  a  very  extensive  study  of  the  history  and 
state  of  the  art  in  the  establishment  of  sandwich  construction  design 
criteria  through  1965,  indicates  that  the  various  stress,  deflection, 
and  stability  analyses  of  sandwich  beams,  plates,  and  shells  are  of  a 
discrete,  rather  than  a  unified,  nature  and  are  somewhat  limited  in 
scope.  For  example,  no  single  analysis  is  referred  to  which  is 


1 


sufficiently  general  to  predict  simultaneously  both  general  instability 
and  face  wrinkling  of  sandwich  plates  or  shells. 

Since  the  compilation  of  the  material  presented  in  Reference  1,  unified 
studies  have  been  completed.  References  2  and  3>  which  are  independent 
and  different  approaches  to  the  same  basic  problem,  deal  with  sandwich 
plates  in  which  both  transverse  shear  and  normal  strain  effects  are  in¬ 
cluded.  In  the  former,  a  variational  approach  is  utilized  and  both  gov¬ 
erning  equations  and  consistent  boundary  conditions  for  prescribing 
either  forces  or  displacements  are  presented.  Applications  are  made 
to  the  stability  in  axial  compression  of  two  types  of  edge-supported 
sandwich  plates  in  References  2  and  4.  In  Reference  3,  governing 
equations  are  presented  for  sandwich  plates  with  dissimilar  faces 
based  on  equilibrium  and  continuity  conditions  developed  in  Refer¬ 
ence  5;  force  boundary  conditions,  established  in  Reference  5  through 
virtual  work  considerations,  are  listed.  An  application  to  buckling 
in  axial  compression  and  bending  of  sandwich  columns  with  dissimilar 
faces  results  in  obviously  incorrect  design  curves.  In  addition,  for 
application  to  plates  (all  edges  supported),  the  buckling  equations 
are  not  consistent  with  the  number  of  boundary  conditions  prescribed. 
Wempner  et  al.,  in  References  5  and  6,  treat  the  problem  of  sandwich 
shells,  the  first  being  a  general  theory  and  the  second  a  specializa¬ 
tion  of  the  general  theory  to  the  specific  problem  of  a  thin  cylindrical 
shell  undergoing  moderately  large  deflections.  However,  when  the  buck¬ 
ling  of  shells  is  treated,  assumptions  are  invoked  in  References  5  and 
6  which  render  the  buckling  equations  unusable  for  studying  all  the 
wrinkling  modes  of  instability.  Also,  in  Reference  5>  the  core- 
constitutive  equations  are  unduly  complex  relative  to  the  core  model 
employed.  In  connection  with  the  work  of  Wempner  et  al.,  it  must  be 
remarked  that,  despite  earlier  presentation,  its  very  late  availability 
in  the  general  literature  (1965)  made  the  theoretical  development  of 
Reference  5  (the  foundation  of  Reference  3)  unknown  to  the  authors 
during  the  development  of  the  present  theory. 
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The  present  study  has  been  undertaken  (1)  as  a  logical  extension  to 
curved  plates  and  shells  of  the  unified  theory  of  Reference  2,  (2)  to 
present  a  set  of  governing  equations,  together  with  consistent  force 
or  displacement  boundary  conditions,  for  the  buckling  and  bending  of 
curved  sandwich  plates  and  shells,  and  (3)  to  correlate  the  theoretical 
predictions  for  the  buckling  of  honeycomb-core  sandwich  shells  in  axial 
compression  with  test  data  available  in  the  literature. 

In  connection  with  item  (1),  the  theory  of  Reference  2  is  obtained  as 
a  special  case  when  both  initial-curvature  and  moderately  large-rotation 
effects  are  omitted. 

Relative  to  item  (2),  the  development  and  form  of  the  governing  equations 
and  boundary  conditions  of  the  present  theory  are  much  less  complex  than 
those  of  References  3  and  5  due  to  the  notation,  structural  models,  and 
method  utilized. 

With  regard  to  item  (3),  there  is  strong  evidence  to  suggest  that  the 
axial  compression  buckling  load  of  a  sandwich  shell  of  practical  con¬ 
struction  and  dimensions,  unlike  the  homogeneous,  isotropic  shell,  is 
given  in  good  approximation  by  the  linear  theory.  The  discrepancies 
that  occur  may  be  as  much  a  function  of  the  changes  in  honeycomb-core 
properties  present  at  the  instant  of  buckling  as  they  are  a  function 
of  initial  imperfections  in  the  geometry  of  the  shell.  Thus,  a  more 
complete  investigation  of  the  problem  suggests  consideration  of  the 
discrete  behavior  of  the  honeycomb  cells  rather  than  the  current  pro¬ 
cedure  of  treating  the  honeycomb  core  as  a  continuum. 


THEORY 


GENERAL  CONSIDERATIONS 

The  basic  approach  taken  in  this  study  is  to  establish  equilibrium 

equations  and  consistent  boundary  conditions  with  the  aid  of  vari¬ 

ational  procedures.  The  equilibrium  and  boundary  condition  equations 
for  the  core  and  faces  (Figure  1)  are  derived  variationally  from  the 
condition  for  a  stationary  character  of  the  separate  total  potentials 
for  the  face  layers  and  the  core  layer.  The  equilibrium  equations  for 

the  core,  applicable  in  the  core  volume,  can  be  integrated  through  the 

thickness  of  the  core.  Then,  by  enforcing  continuity  of  displacements 
at  the  interfaces  between  the  core  layer  and  the  face  layers,  a  final 
set  of  governing  equations  and  consistent  boundary  conditions  is  ob¬ 
tained  in  only  two  surface  coordinates,  comparable  to  the  usual  two- 
dimensional  formulation  associated  with  thin  shell  theory. 

FACE  CONSIDERATIONS 

The  face  layers  are  treated  as  thin,  isotropic  shells  on  the  basis  of 
the  Kirchhoff-Love  assumptions.  Rotations  about  the  normal  to  the 
middle  surface  are  neglected  with  respect  to  rotations  about  direc¬ 
tions  tangent  to  the  middle  surface. 


The  middle  surface  strain  displacement  relations  and  the  curvature 
changes  are  (Reference  7) 
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As  observed  in  Appendix  I,  the  strain  energy  of  a  linearly  elastic  iso¬ 
tropic  medium  for  which  the  Kirchhof f-Love  approximations  are  valid  is 
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and  i  takes  on  the  values  1  and  2  for  the  upper  and  lower  faces, 
respectively. 
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the  potential  of  applied  loads  for  one  of  the  faces  is 
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where  the  upper  sign  applies  to  face  1  (i  =  1),  and  the  lower  sign 
applies  to  face  2  (i  =  2). 

Application  of  the  minimum  total  potential  energy  principle  requires 
that 
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The  actual  variation  is  carried  out  in  Appendix  I.  The  resulting  Euler 
equations  of  equilibrium  are 
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The  boundary  conditions  consistent  with  the  equilibrium  equations  are 
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where  c^  ,  c  ,  ,  d^2  ,  etc.,  are  transformation  symbols  utilized 

to  express  the  boundary  conditions  in  terms  of  generalized  displacement 
components  corresponding  to  a  curvilinear  coordinate  system  along  and 
normal  to  the  edge  of  the  shell. 


CORE  CONSIDERATIONS 

The  core  description  utilized  herein  is  representative  of  that  associated 
with  honeycomb-sandwich  construction.  For  such  a  core,  it  is  reasonable 
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to  assume  that: 

1.  The  inplanc  shear  and  extensional  stiffnesses  are  negligible 
with  respect  to  transverse  stiffness  quantities,  corresponding 
to  a  state  of  stress  in  which  rr^  =  rr0  =  0 

?.  In  view  of  item  1,  it  is  justifiable  to  linearize  the  rele¬ 
vant  strain-displacement  relations  for  the  core. 

3.  The  transverse  stiffnesses  per  cell  are  constant  through  the 
thickness. 

The  first  assumption  is  obvious  in  the  case  of  a  ree  core  of  cellular 
structure.  When  the  core  is  supported  by  faceplates,  the  assumption 
is  reasonable  for  thick  cores  (c  >  a)  and  for  laterally  loaded  sand¬ 
wich  panels  with  relatively  low  lnplane  loads.  However,  for  very  thin 
panels  (for  example,  c  <  ^  )  loaded  in  the  plane  of  the  structure, 
the  assumption  may  need  reconsideration,  especially  when  face-wrinkling 
phenomena  are  being  studied. 

An  element  of  the  sandwich  shell  under  consideration  is  shown  in  Fig¬ 
ure  1.  Tht  middle  surface  of  the  undeformed  core  (£  =  0)  is  utilized 
as  the  reference  surface  of  the  shell.  Orthogonal  curvilinear  coor¬ 
dinates  and  are  taken  along  the  curvature  directions.  The 

principal  radii  of  curvature  are  positive  as  shown  in  the  figure.  The 
coordinate  f  is  measured  normal  to  the  surface. 

On  the  basis  of  the  core  model,  the  relevant  constitutive  equations  for 
the  resulting  antiplane  state  of  stress  are 
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In  Equation  (8),  the  stiffness  functions  for  a  curved  element  are  ap¬ 
propriately  designated  as 
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where  a2  ,  a2  ,  f  ,  ,  and  C2  are  the  values  of  these  quantities 

when  £  =  0  . 

The  linearized  strain  displacement  relations  are 
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With  the  use  of  Equations  (6)  and  (10),  it  is  shown  in  Appendix  II  that 
the  strain  energy  of  the  core  is 
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For  the  edge  load  q  and  the  interface  load  system  with  components 
q^  ,  q2  ,  and  q^  (Figure  2),  the  potential  of  applied  loads  is 
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When  the  minimum  total  potential  energy  principle  is  applied  to  the 
core  in  the  form 
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the  equilibrium  equations  and  b>.  mdary  conditions  are  obtained  as  shown 
in  Appendix  II.  The  equilibrium  equations  are 
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The  boundary  conditions  consistent  with  these  equilibrium  equations  are 
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INTERFACE  CONTINUITY  CONSIDERATIONS 


At  the  present  point,  the  governing  equations  and  boundary  conditions 
have  been  derived  for  the  core  and  faces  considered  as  separate  media. 
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In  order  to  ensure  continuity  of  deformation  in  the  composite  sandwich 
shell,  the  components  of  the  core  and  face  displacements  must  satisfy 
the  conditions 
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As  demonstrated  in  Appendix  II,  the  satisfaction  of  conditions  (l6a,b, 
and  c)  subsequent  to  integration  of  the  core  strain-displacement  equa¬ 
tions  gives  rise  to  the  two  interface  compatibility  equations 


o  (17) 


15 


and 


1  \  2p0/  M  ?p0/  «0  ?  6P_ 


1  C  +  t  C  +  3t?C 

+  —  w  - -  +  - —  -  ye 

a2  ?’*2  ?  6p? 


c  G1  -  GP  - 

— —  (<Vi)  t  +  rr  (<v9) 

12  Fa2  |aja?  ,51 


As  is  also  shown  in  Appendix  II,  additional  results  stemming  from  the 
integration  of  the  core  equations  permit  the  interface  loads  q^  ,  q2  , 
and  q^  and  the  boundary  conditions  for  the  edge  of  the  core  to  be 
expressed  in  terms  of  the  face  displacement  components  u^  ,  v^  ,  an(* 
w^  and  the  core  shear  angles  y ^  and  y 2  . 

GOVERNING  EQUATIONS  AND  BOUNDARY  CONDITIONS 

With  terms  usually  eliminated  from  shallow-shell  theory  denoted  by  an 
asterisk,  the  final  set  of  governing  equations  describing  the  behavior 
of  the  sandwich  shell  is 
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The  boundary  conditions  in  final  form,  which  permit  either  generalized 
forces  or  displacements  to  be  prescribed,  are 
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The  governing  equations  and  boundary  conditions  given  by  Equations  (17), 
(18),  (19),  and  (20)  are  sufficient  to  determine  the  bending,  buckling, 
and  postbuckling  behavior  of  sandwich  shells  under  prescribed  loads  and 
boundary  conditions. 


Special  forms  and  an  application  of  these  equations  are  treated  in  sec 
tlons  to  follow. 
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Shallow  Shell  Equations 

In  many  instances,  the  contribution  of  inplane  displacements  to  the 
rotations  m  and  o>2  of  the  face  layer  is  small.  Especially  in 
the  case  of  shallow-sandwich  shells  and  shallow  regions  of  nonshallow 
shells,  the  expressions  for  the  rotations  can  be  simplified  to  give 
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If  these  rotation  expressions  are  utilized,  then  terms  marked  by  an 
asterisk  have  to  be  omitted  from  the  governing  equations  and  boundary 
conditions.  In  such  a  case,  the  equations  for  the  faceplates  are  es¬ 
sentially  those  commonly  referred  to  as  the  von  Karman-Donnell  type. 


Thin  Shell  Equations 

In  cases  where  the  total  shell  thickness  (c  +  t^  +  t^)  is  negligible 
with  respect  to  each  of  the  principal  radii  of  curvature  (p^,p^)  , 

terms  of  the  order  (c  +  t^  +  tg)/p  or,  in  view  of  the  Kirchhof f-Love 
assumption  (t/p)  «  1  ,  of  the  order  c/p  can  be  neglected  with  re¬ 

spect  to  unity.  As  a  further  consequence,  the  differences  in  the  metric 
components  (Q^jOig)  and  the  radii  of  curvature  (Pj,Pg)  of  the  two 
face  layers  must  be  neglected.  With  these  simplifications,  the  gov¬ 
erning  shallow-shell  equations  and  boundary  conditions  of  the  previous 
section  reduce  to 
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These  equations  for  the  sandwich  configuration  in  which  F  -♦<»  are 
those  of  Reference  8.  Further,  if  the  bending  stiffness  of  the  faces 
is  neglected  (membrane  faces),  then  the  equations  represent  the  gen¬ 
eral  form  of  the  classical  theory  for  sandwich  shells  of  Reference  Q. 
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APPLICATION  TO  COMPRESSIVE  BUCKLING  OF  THIN  CIRCULAR  CYLINDERS 


GOVERNING  EQUATIONS  AND  BOUNDARY  CONSIDERATIONS 

The  present  theory  is  applied  to  the  problem  of  determining  critical 
loads  for  a  uniformly  axially  compressed  cylinder  of  honeycomb  sand¬ 
wich  construction.  Attention  is  directed  especially  to  establishing 
the  effect  of  core  stiffness  on  the  distinct  instability  modes  and  the 
interactions  of  these  phenomena.  The  cylinder  is  shown  in  Figure 


In  Appendix  III,  the  shallow-shell  equations  are  specialized  for  the 
case  of  cylindrical  coordinates.  Uniform,  radially  symmetric  deforma¬ 
tions  are  assumed  to  prevail  under  the  axisymmetrlc  load  up  to  a  criti¬ 
cal  value;  that  is,  edge  restraints  are  considered  not  to  influence  the 
deformation  prior  to  buckling.  To  determine  the  possible  bifurcation 
points  of  the  load-shortening  relationship,  infinitesimal  deformations 
measured  from  the  unbuckled,  compressed  state  are  used.  The  special¬ 
ized  governing  equations  (Equations  (119)  of  Appendix  III),  in  terms 
of  these  deformations,  are 
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The  appropriate  boundary  conditions  for  the  specialized  governing  equa¬ 
tions  (Equation  (119)  of  Appendix  III)  are 
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As  pointed  out  in  the  development  of  Appendix  III,  the  thickness  of  the 
cylinder  has  been  considered  negligibly  small  with  respect  to  the  radius; 
thus,  the  compressive  load  is  distributed  evenly  over  the  two  faces. 

In  addition  to  the  general  instability  type  of  buckling,  the  phenomenon 
of  very  short  wave  deformations  (face  wrinkling)  receives  special  atten¬ 
tion  herein.  The  length  of  a  cylinder  is  supposed  to  be  large  enough  as 
compared  to  the  wavelength  of  any  typical  deformation  pattern  to  allow 
the  effect  of  boundary  conditions  on  the  magnitude  of  the  critical  loads 
to  be  ignored.  The  set  of  nine  boundary  conditions  (Equation  (25))  gives 
rise  to  a  very  large  number  of  possible  combinations  of  generalized  forces 
and  displacements.  However,  although  boundary  conditions  are  of  decisive 
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importance  for  edge  deformation,  their  importance  decreases  as  the  num¬ 
ber  of  repeated  deformation  patterns  along  the  length  of  the  cylinder 
increases. 

EQUIVALENT  SINGLE  GOVERNING  EQUATION 

It  is  convenient  to  consider  a  single  partial-differential  equation  in 
one  of  the  unknowns  rather  than  the  eight  simultaneous  equations  in 
eight  unknowns.  In  this  way,  the  relative  importance  of  the  stiffness 
parameters  and  geometric  quantities  can  be  assessed,  and  the  character¬ 
istic  equation  for  the  "classical"  solution  can  be  extracted  directly; 
in  addition,  the  specialization  to  previous,  but  less  general,  theories 
for  plates  and  shells  can  be  readily  accomplished. 

The  governing  equations  can  be  written  in  the  form 

A  r  =  0  (26) 

where  A  is  the  matrix  of  linear  operators  and  r  is  a  column  vector 
whose  components  are  the  functions  ,  u^  ,  v^  ,  v?  ,  w^  ,  w2  ,  cy^  ,  and  cy^  . 

The  matrix  A  is  diagonalized  when  postmultiplied  by  its  adjoint  to 
yield 

A  IjA^ir  =  |A|  I  (27) 

where  ||A^j||'  is  the  adjoint  of  A  or  the  transposed  matrix  of  the 
cofactors  of  the  elements  of  A  ,  and  I  is  the  unit  matrix.  | A | 
is  the  determinant  of  A  .  Instead  of  Equation  (26),  the  equation 

|A |  I  r  =  0  (28) 

may  be  solved;  that  is,  a  nontrivial  solution  to  the  equation,  say,  for 
the  radial  displacement  w^  , 

|A|  Wj  =  L(wx)  =  0  (29) 

must  be  found.  The  symmetric  operator  matrix  A  is 


29 


1 


30 


The  linear  operator  L 
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The  order  of  the  highest  derivative  in  this  operator  is  18;  the  number 
of  independent  conditions  that  can  be  imposed  at  the  boundary,  in  terms 
of  generalized  forces,  displacements,  or  combinations  thereof,  is  9  • 
Equation  (31)  relative  to  existing  equations  of  thin  plates  and  shell 
theory  is  discussed  in  a  subsequent  section. 

It  may  be  noted  that  the  operator  in  the  governing  equation  (Equation 
(29) )  can  be  written  as 


L  =  L0  +  (if  L1  w 

As  terms  of  the  order  c/R  have  been  neglected  with  respect  to  unity, 
the  second  term  on  the  right-hand  side  of  Equation  (32)  can  be  ignored 
if 
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'.n  estimate  of  the  relative  importance  of  the  last  part  of  the  operator 
L  is  made  in  conjunction  with  the  analyses  of  the  axisymmetric  and  non- 
axisymmetric  buckling  modes  carried  out  in  the  following  sections. 

Axisymmetric  Buckling 

For  the  case  of  axisynmetric  behavior,  the  deformations  v^  ,  v^  ,  and 
y g  ,  as  well  as  all  derivatives  with  respect  to  cp  ,  vanish  and  the 
resulting  governing  operator  becomes 
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Thus,  for  the  case  of  axlsymmetric  buckling,  the  18th  order  governing 
equation  (Equation  (29))  can  eplaced  by  a  14th  order  equation.  On 
the  basis  of  the  same  considerations,  the  number  of  boundary  conditions 
is  reduced  from  9  to  7  .  Specifically,  the  boundary  conditions  given 
by  Equations  (25c)  and  (25d)  become  trivial. 


A  solution  to  the  governing  equation  L  (r)  =  0  ,  corresponding  to 

3 

"classical"  simple  support  of  the  faceplates  (vanishing  radial  edge 
displacement  and  edge  moment),  is 
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The  forms  of  u^  ,  u0  ,  and  > ^  are  implied  by  the  relationships  that 
exist  between  these  quantities  and  the  radial  displacements  w^  and  w. 
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in  the  original  governing  equations.  The  "implied"  boundary  conditions 
for  the  case  of  "classical"  simple  support  of  the  faces 
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With  the  use  of  the  dimensionless  parameters 
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a  characteristic  equation  results  as 
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Equation  (36)  is  of  the  form 


Ak  -  Bk  +  C  =  0 


(37) 


where  A  ,  B  ,  and  C  are  polynomials  in  the  wavelength  parameter  p  . 

2 

As  A  ,  B  ,  C  ,  and  B  -  4AC  are  positive  definite  quantities,  two 
real-valued  positive  roots  are  obtained  for  each  value  of  the  wavelength 
corresponding  to  different  values  of  tbj  ratio  W^/W^  '  relat*-ve 

importance  of  the  last  term  of  the  left-hand  side  of  Equation  (36)  can 
be  readily  determined  by  examining  the  coefficients  B  and  C  of 
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Equation  (57).  With  N  defined  as  fc  v  / ( 1 - v )  and  irrelevant  terms 
omitted,  the  expressions  become 
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is  of  the  order  of  magnitude  c/R  for  typically  large  values  of  E/2F 

(for  example,  2  X  10  )  and  t/R  (for  example,  5  X  10  ^  )  and  is 

therefore  to  be  neglected  with  respect  to  unity.  Hence,  the  con- 

O  p 

tribution  of  N  to  the  last  term  in  both  B(3  )  and  C(f3  )  is 

negligible.  The  effect  of  N  in  the  coefficient  [3  4  in  C([32) 


is  certainly  small,  as 

N 

cpv 


i?  ‘  “°(i) 


These  and  other  simplifications  of  the  same  nature,  although  not  essen¬ 
tial  for  a  numerical  computation  of  critical  values  of  the  load  and  the 
wavelength,  allow  an  analytical  evaluation  of  both  in  closed  form.  When 
the  N-term  is  negligible,  the  characteristic  equation  simplifies  to 


(32k  -  (Pk  +  rp  +  £  j 


16^  4  22) 

,  Y  +  0pp  +  t  &  j  k 


-  j  3(1 +  c)2  04  +  ^-  +9p32  +  e2j(p4+fJ)J 


0  (38) 


The  characteristic  equation  (Equation  (36))  and  its  simplified  form 
(Equation  (38))  are  solved  for  the  load  parameter  k  by  using  typical 
values  of  the  curvature,  relative  thickness,  and  stiffness  parameters. 
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The  load  parameter  k  versus  the  wavelength  parameter  p  ,  considered 
to  be  a  continuous  variable,  is  plotted  in  Figure  4.  In  Figure  5,  a 
similar  solution  is  presented  for  values  of  the  various  parameters 
corresponding  to  low  values  of  F  and  to  assess  qualitatively 

the  buckling  of  a  relatively  weak-core  sandwich  shell.  For  compari¬ 
son,  corresponding  solutions  based  on  existing  theories  (References 
11  and  8)  are  shown  also.  No  quantitative  distinction  can  be  made 
between  the  solutions  of  the  complete  and  simplified  quadratic  charac¬ 
teristic  equations;  thus,  only  a  single  curve  is  evidenced  for  each 
branch  of  the  solution. 

The  characteristic  equation  is  of  the  form 


(k  -  k1)  (k  -  k2)  =  0 


(39) 


where 


?  P  fp  +  i 
k  .(r)  =  P*  +  — =- 
P 


(*0) 


and 


kg(P  )  = 


3(1  +  e)2  p2 


+  0pP2  +  e2 

cp 


i 

+  p2  +  — 

p2 


(41) 


Q 

Stationary  values  of  k^  or  k^  occur  whenever  Sk/8(p  )  =  0  ,  and 

a  critical  situation  is  reached  when  the  stationary  value  of  k  is  a 

2  2  2  p 

minimum  (P  k/4(p  )  >  0)  .  As  p  is  a  real-valued  quantity,  p 

is  always  nonnegative  and,  in  practice,  always  positive.  The  root 
kj  has  only  one  stationary  value;  a  minimum  occurs  when 


M, 


=  (v  +  O 


1/4  1/4 

'  *  cp  ' 


(42) 
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and  the  corresponding  critical  load  is  given  by 


(kcr )  =  2  \.cp  +  £  ~  2  x/cp"  (43) 

since,  as  mentioned  earlier,  f^/cp  may  be  neglected  in  comparison  with 
unity.  The  root  usually  has  more  than  one  stationary  value  in 

the  p  interval  under  consideration.  This  may  be  seen  in  Figures  4 
and  5  •  For  the  more  practical  set  of  core  parameters,  minimums  occur, 

corresponding  to  two  very  distinct  values  of  p 


The  longwave  instability  mode  (p  «  l)  can  be  analyzed  by  neglecting 
terms  which  are  obviously  small  due  to  the  small  size  of  6  .  Equa¬ 

tion  (4l)  then  simplifies  to 


k 


2 


3(1  +  ef  p2  £ 
0PP2  +  c2  p2 


(44) 


From  this  expression,  a  critical  wavelength  parameter 


c(l  +  c) 


\/3£  -  0P^ 


-1/2 


(^5) 


and  the  corresponding  critical  load  parameter 


(kcr)  -  2(1  +  0  y  /*  ~  OP *J<2  (46) 

cl 

are  determined.  It  may  be  noted  that  the  terms  neglected  in  Equation 
(44)  are  those  originating  from  the  inclusion  of  the  effects  of  face 
bending  stiffness  and  the  core  normal  strain.  Obviously,  the  result 
obtained  corresponds  to  the  classical  solution  of  Reference  11;  there¬ 
fore,  the  critical  wavelength  and  load  have  been  assigned  a  subscript 
"cl"  for  this  case. 
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For  the  shortwave  instability  mode,  g  is  of  the  order  of  magnitude  of 
unity  or  even  larger.  Equation  (kl)  may  then  be  simplified  to 


3(1  +  c)2 

k  ~  — -  +  P 

p/cp  +  0  P 


The  critical  wavelength  and  associated  load  parameter  are 


(1  +  0 


3cp  -  0pq' 


1/2 


(47) 


(48) 


and 


( k  )  =  2(!  +  c)  3q^  -  0Pcp  (49) 

2 

The  results  presented  in  this  section  are  discussed  in  detail  in  the 
subsequent  section  entitled  "RESULTS  AND  DISCUSSION". 

Nonaxisymmetric  Buckling 

For  the  case  of  nonaxisymmetric  buckling,  a  solution  of  the  general 
equation  (Equation  (28))  is  sought.  Such  a  solution,  corresponding 
to  "classical"  simple  support  of  the  faceplates  (vanishing  edge  dis¬ 
placement  and  edge  bending  moment),  is  available  in  the  form 


U1 " 

U^  cos  7Tx/X  sin  nrp 

U2 

Ug  cos  7Tx/X  sin  ncP 

V1 

V  sin  7 tx/X  cos  ncp 

V2 

V  sin  ttx/X  cos  nqi 

W1 

W  sin  7Tx/X  sin  nrp 

W2 

sin  7Tx/X  sin  ncp 

c0 

cl  cos  7Tx/x  sin  ntp 

_C ;  2  _ 

cF^  sin  7rx/X  cos  nrp 

31 


(50) 


} 


The  forms  of  ,  7^  ,  and  7 ^  are  implied  by  relation¬ 

ships  that  exist  between  these  quantities  and  the  radial  displacements 
w^  and  w  in  the  original  governing  equations.  The  "implied"  boundary 
conditions  for  the  case  of  "classical"  simple  support  of  the  faceplates 


are 


=  D 


W2  =  ° 


w,  +  vw, 

1 J  XX  l>cpcp 


W-  +  vw, 

2, xx  1 ,  qxp 


=  0 


-  0 


(51a) 

(51b) 


w. 


U1,X  +  V 


=  0 


(51c) 


V!  =  0 


(51d) 


w. 


U2,x  +  V 


=  0 


(51e) 


v2  -  0 


(5  If ) 


c 

6f 


"1  +  “2  +  ~  Vi,x  “  0 


(51s) 


With  the  use  of  dimensionless  parameters  previously  defined  and  the 
buckle  aspect  ratio  parameter 


P  = 


2^2 
n  X 

o  p 

7T  R  \  X 


r 


a  characteristic  equation  is  obtained  as 
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p 
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(p  +  2)  |  +  —  P  (l  +  PP  )! 


1-V  1+V 


l-v  P 


=  o  (52) 


This  equation  reduces  to  l.._  characteristic  equation  for  the  case  of 
axisymmetric  buckling  (Equation  (36))  when  the  wave  aspect  ratio  p 
is  set  equal  to  zero.  For  nonaxisymmetric  buckling,  p  is  positive 
definite  and  the  last  part  of  the  characteristic  equation,  marked  by 


Ul 


an  asterisk,  can  be  ignored  if  the  condition 


£/cp  ^  0(c/R) 

is  satisfied. 

Then,  Equation  (52)  reduces  to  the  form 

(k  -  k1)  (k  -  k2)  =  0 

where 

o  o  ^  +  Cp 

k  =  p  (l  +  pr  +  — 

p 

and 


(53) 


(54) 


(55) 


k  -  p2(l+P)2  +  -7 - 5 

2  P2(1+p)2 


3(l+e)2  P6(l+p)2 

~~  +  -  P2(1+PP2) 
1-V  P 

0p(l+p)p6 

f.liiL Ze!t£p2(I+p)] 

1-V  l+p  P  v 

+ 

e2P^  +  ^5  (1+p)2 

pp2) 

(56) 


The  first  factor  on  the  left-hand  side  of  Equation  (54)  leads  to  the 
critical  value  of  the  wave  parameter 
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(cp  +  O 


1/4 


1  +  p  ,  1  +  p 

and  the  associated  critical  load  parameter 


Kr). 


(57) 


p=0 


(k  )  =  2(1  +  p)  fj'  +  f  »  2(1  +  p)  fp  -  (l+p) 

c  r  - 


(ker). 


p=0 

(58) 
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Since  p  is  positive  definite,  the  critical  wavelength  and  the  criti¬ 
cal  load  are  larger  than  In  the  ease  of  nxlnyrmotr lc  buckling  (p  0)  . 


The  second  root  of  Equation  (5M>  k  ,  la  rewritten  by  dividing  the 
last  terms  by  the  nonzero  factor 


0  2  2 

—  +  -  fr(l  ♦  pi r) 

l-v  t> 


Then, 


0?(l+p)?  ♦ 
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2  es(u,)  .zl.hili 
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The  factor  f  depends  on  the  ratio  of  the  core  shear  moduli  pc 


(59) 


(60) 


and 

It 


Is  equal  to  unity  for  the  case  of  an  isotropic  core  (p  -  1) 

2 

can  be  observed  that  the  dependence  of  f  on  0  Is  not  large,  since 
p  =  Cg/ Cj  possesses  values  In  the  Interval  0.6  to  1.5  for  practi¬ 

cal  honeycomb-core  materials.  The  wave  aspect  ratio  parameter  p  «*  (X  /X  ) 


Is  less  than  or  about  0.25  for  longwave  deformations. 
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In  particular,  the  ratio 
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Is  ■null  with  respect  to  unity.  If  p  ■  1  (Isotropic  core),  thtn  tho 
ratio  vanlshaa.  For  example,  whan 


_  * 
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For  very  short  wavelengths, 


f 


1  +  p 

- 5  «nd 

1  +  PP* 


-*  0 


p 

Upon  neglect  of  the  dependence  of  f  on  (J  end  with  the  assumption 

o 

that  P  Is  small  compared  to  unity  (longwave  buckling),  the  critical 
values  for  the  wavelength  and  load  parameters  are  obtained  as 


c(l+c)  (1+p) 


(62) 


2(1* )  Vx/e2 


1*9 


(63) 


0pCA 
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for  a  given  value  of  p  .  Equation  (63)  above  that  when  tha  core 
ahear  atlffnasa  approachea  Infinity  (0  -*  0)  ,  then 


-  2(1*)  Vx/c2' 


Thla  raault  corraaponda  to  nonaxlsymatrlc  claaalcal  buckling  of  a 
hoaogeneoue  lao tropic  ahell  and  la  Identical  to  tha  raault  obtained 
for  tha  axlayamtrlc  soda.  However,  when  tha  ahear  atlffneaa  la 
finite,  Equations  (63)  and  (U6)  (axlsymetrlc  buckling  caaa)  Indi¬ 
cate  that  tha  tvo  claaalcal  aandvlch  shall  raaulta  are  different. 
Indeed,  for  practical  values  of  tha  core  shear  stiffness  parameter 
(0.0  <  p  <  1.5)  ,  tha  nonaxlayaanetrlc  buckling  n»de  leads  to  a 
slightly  different  value  of  the  buckling  load.  The  applicability 
of  these  « qua t Ions  for  long  shells  to  shells  of  finite  length  Is 
discussed  later. 


For  a  shortwave  buckling  node  (3  -0(1))  ,  the  critical  wavelength 

par aae ter 


(l*)  -  fepep  x/2 


(6k) 


and  the  corresponding  critical  load  parameter 
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2(1*)  '  3^  ’  f0P«P 
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2(1*)  -.3<p 
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since  for  very  short  wavelengths, 


Up 

f  -*  2 
1+pp^ 

Equation  (65),  when  compared  with  Equation  (49)  for  the  axlsynaetrlc 
case,  Indicates  a  higher  value  of  the  critical  load  parameter  when 

(1+p)2  2(Uc)  5qT 

— T  <  1  ♦ - 

l+pp  0  pq? 

In  most  cases  this  condition  will  be  met;  a  more  complete  discussion 
follows  In  the  section  entitled  "RESULTS  AND  DISCUSSION". 


(66) 
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RESULTS  AND  DISCUSSION 


THEORY 

The  general  theory  developed  In  thli  Investigation  la  given  by  Equations 
(17),  (16),  (19),  and  (20).  These  relations  represent  the  governing 
equations  and  associated  boundary  conditions  which  can  be  used  to  de¬ 
scribe  the  bending,  buckling,  and  postbuckllng  behavior  of  sandwich 
shells  (and  curved  plates).  This  theory  Is  applicable  to  srndwlch 
shells  which  arc  constructed  of  thin,  Isotropic  face  layers  and  a 
honeycomb- type  core  layer;  It  Is  considered  to  be  the  most  complete 
and  consistent  theory  developed  to  date  for  such  a  configuration. 

Of  the  many  sandwich  shell  theories  which  exist  lr  the  literature, 
the  recent  studies  (References  5  and  6)  by  Wenpner  et  al.  are  the 
moat  advanced  studies  with  which  the  present  development  can  be  com¬ 
pared.  It  would  appear  from  a  cursory  examination  that  the  two  the¬ 
ories  essentially  describe  the  same  problem;  nevertheless,  It  la 
evident,  as  will  be  disclosed,  that  there  are  narked  and  Important 
differences  relative  to  both  development  and  applicability. 

A  direct  comparison  with  the  work  of  Wempner  et  al.  reveals  that: 

1.  The  treatment  of  the  face  layers  Is  equivalent. 

2.  The  core  considerations  utilized  In  Reference  5  are  baaed 
on  a  model  In  which  a  contradiction  exists  between  the 
linearized  equilibrium  equation  employed  and  the  nonlinear 
strain-displacement  relation  assumed.  That  is,  the  core 
Is  assumed  to  be  in  a  linear  state  of  antiplane  stress 
while  the  kinematic  relations  tacitly  Imply  a  more  com¬ 
plex  behavior. 

3.  Even  though  both  theories  take  into  account  the  effect  of 
core  thickness  with  respect  to  radii  of  curvature,  the  the¬ 
ory  of  Reference  3  does  not  account  for  the  fact  that  a 
honeycomb  core  (which  is  within  the  weak  core  family  for 
which  the  theory  is  stated  to  be  valid)  does  not  possess 
the  same  core  shell  area  at  the  the  two  intertaces.  For 
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Chin  sandwich  sheila  (c/p  negligible),  this  consideration  Is 
noC  Important.  However,  for  thlck-corc  sandwich  shells,  the 
only  example  which  la  Illustrated  In  Reference  5,  the  con* 
eldoratlon  la  Important  and  rendera  the  results  questionably 
applicable  to  a  honeycomb-core  sandwich  or,  In  fact,  to  any 
thick-core  sandwich  In  which  an  Initially  flat  core  medium 
la  fonaed  Into  a  curved  aurface. 

4.  In  Reference  6,  which  specifically  treats  thin  aandwlch  shells, 
the  final  equations  presented  are  Inapplicable  to  the  complete 
study  of  the  buckling  of  such  shells,  since  a  gross  assumption 
la  Invoked  that  precludes  the  occurrence  of  the  antlsynnetrlc 
face-wrlnkllng  mode  of  Instability. 

5.  The  boundary  conditions  presented  by  Wempner  et  al.  are  de¬ 
veloped  Independently  of  the  governing  equations  and  are 
stipulated  In  terns  of  edge  forces  only.  The  present  theory 
develops  both  governing  equations  and  force-displacement 
boundary  conditions  within  the  consistent  framework  of  a 
complete  variational  treatment. 

The  governing  equations  (Equations  (IT),  (18),  and  (19))  of  the  present 
theory  can  be  shown  to  be  equivalent  to  a  single  l8th-order  partial 
differential  equation  In  any  one  of  the  variables  (for  example,  the 
normal  displacement  of  a  face  layer).  The  associated  boundary  con¬ 
ditions,  given  by  Equation  (20),  are  consistent  with  the  l8th-order 
system  since  9  boundary  conditions  on  either  generalized  force  or 
displacement  can  be  prescribed  at  an  edge  of  the  shell.  A  special 
case  of  the  general  equations  Is  applicable  to  the  buckling  of  a  thin, 
circular  cylindrical  shell  (see  Equation  (31)).  This  equation  has  been 
obtained  by  systematically  reducing  the  eight  governing  equations  for 
the  thin,  circular  cylindrical  sandwich  shell  in  axial  compression 
(Equation  (24))  through  the  use  of  linear  operators. 

Before  undertaking  the  discussion  of  the  results  of  the  application  of 
the  present  theory  to  the  compressive  buckling  of  thin,  circular  cylin¬ 
drical  shells,  it  is  of  interest  to  show  the  relationship  of  the  present 


single-governing  equation  of  the  18th  order  with  previoua  equal Ions 
utilised  to  study  the  buckling  of  flat  plates  and  circularly  curved 
plates  and  cylinders,  Including  the  associated  requirement  for  suf¬ 
ficient  boundary  conditions  to  effect  complete  solutions.  With 
reference  to  Table  I,  It  can  be  seen  that  the  present  l8th-order 
equation,  requiring  9  boundary  conditions  per  edge,  reduces  to  the 
7  other  cases  Indicated  when  the  prescribed  simplifying  assumptions 
ere  enforced.  From  the  table,  It  is  noted  dramatically  that  the  trend 
toward  higher  performance  structures  evidences  increasing  complexity 
with  regard  to  theoretical  considerations.  This  point,  although  the 
result  of  theoretical  considerations,  cannot  be  casually  dismissed  by 
Individuals  concerned  with  the  design,  development,  analysis,  and 
testing  of  aerospace  components  which  reflect  a  strong  structures- 
materlals  Interface  aimed  at  achieving  Increased  structural  efficiency. 

AXISYMMETRIC  BUCKLING  OF  THIN.  CIRCULAR  CYLINDRICAL  SANDWICH  SHELLS 

The  solution  of  Equation  (29)  for  the  case  of  axlsymmetrlc  buckling 
leads  to  the  buckling  criterion,  Equation  (36),  which  Is  presented 
graphically  In  Figure  4.  The  figure  shows  the  variation  of  the 
buckling  load  parameter  k  with  wavelength  for  the  cylinder  de¬ 
scribed  by  the  parameter  values  noted.  As  mentioned  earlier,  these 
parameter  values  reflect  an  actual  cylinder,  the  test  results  for 
which  are  described  In  Reference  10.  Figure  4  reveals  that  three 
distinct  mlnlmums  occur  for  practical  values  of  the  core  stiffness 
parameters.  Each  minimum  corresponds  to  a  different  buckling  mode; 
namely,  one  relatively  long-wave  antisymmetric  mode  and  two  relatively 
short-wave  modes  corresponding  to  antisymmetric  and  symmetric  buckling. 
When  the  case  of  a  shell  with  very  low  values  of  the  core  moduli  is 
considered,  a  different  behavior,  leading  to  consonance  of  the  two 
usual  anti8yimetric  modes,  is  evidenced  (Figure  5)*  Although,  on  the 
basis  of  weight  considerations,  such  a  structure  might  appear  to  be 
advantageous,  the  low  shear  stiffness  and  the  correspondingly  large 
shear  deformations  would  render  the  composite  impractical. 


TABLE  I 

COMPARISON  OF  PLATE  AND  SHELL  THEORIES  WITH  RESPECT  TO  ORDER  OF  GOVERN IN 
PARTIAL  DIFFERENTIAL  EQUATION  AND  BOUNDARY  CONDITIONS  REQUIRED  PER  EDGE 
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Figure  U.  Relation  Between  Load  Parameter  k  and  Wavelength 
Parameter  0  for  a  Sandwich  Cylinder  of  Practical 
Properties. 
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Figure  5.  Relationship  Between  Load  Parameter  k  and  Wavelength 
Parameter  3  for  a  Hypothetical  Weak-Core  Sandwich 
Cylinder. 
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The  character  of  the  different  buckling  modes  (Figure  •»)  can  be  readily 
recognized  If  a  deformation  vector 


u 

U1  +  u? 

A* 

u 

U1  -  U? 

A* 

w 

1 

e  — 

2 

w,  +  w. 

1  2 

w 

*1  -  Wg 

C7 

_?C>1 

(67) 


la  considered.  By  performance  of  elementary  row  operations  on  the 
matrix  A  (Equation  (26)),  and  after  substitution  of  the  solution 
given  by  Equation  (3U),  this  vector  can  be  written  as 


ll  p 

(3  -  kfT  +  (  +  qp 


-JV/2 


(1-v  ) 


(1-K 


6(l-v?)1/2  0  /  u  2 

-  - —  p4  -  V8T  +  f  + 

7TV  c  y  ' 


*) 


!  — ( 

7T  P(l+e)  V 


p1*  -  kp2  +  f 


Op 


'  ( 1+0  3r 


) 


—  (p1*  -  kp2  +  r  )  (p14  -  kpg  +  C  +  q>) 


(66) 
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For  a  solution  to  Equation  (4o),  only  the  components  u  and  w  arc 
nonzero.  Since  the  ratio 

U  7TV  C  f  VC 

I  — - X-  -  *  ~  *  ~  (69) 

W  12(l-v  j  3  7jp  R 

and,  usually,  since  3  >  1  ,  the  relative  face  displacement 

is  seen  to  be  the  significant  deformation.  Thus,  the  buckling  mode 
is  characterized  by  essentially  symmetric  displacements  of  the  faces. 

The  wavelength  of  this  instability  mode  is  of  the  order  of,  or  smaller 
than,  the  core  thickness  (Equation  (1+2));  this  mode  is  referred  to  as 
the  symmetric  wrinkling  mode. 

On  the  basis  of  a  comparison  similar  to  that  made  in  Equation  (69), 
the  solutions  corresponding  to  Equations  (44)  and  (47)  are  recognized 
as  being  essentially  antisymmetric  in  nature  (w  »  w  ,  u  »  u) 

These  characteristics  are  more  pronounced  in  the  case  of  the  longwave 
buckling  mode  than  in  the  wrinkling  mode.  The  former  has  already  been 
related  to  the  classical  instability  mode  (Equation  (46));  the  latter 
is  referred  to  as  the  antisymmetric  wrinkling  mode. 

In  Figures  4  and  5>  the  curves  for  the  antisymmetric  solution  based  on 
the  classical  theory  of  Reference  9  approaches  an  asymptotic  value  for 
large  values  of  the  wavelength  parameter.  Equation  (44),  which  is  valid 
also  if  F  — )  00  ((p  — >  00  )  ,  reveals  that  for  increasing  (3  ,  the  load 
parameter  k^  approaches  the  classical  theory  value  3/0P  if  e  =  0 
(membrane  faces).  The  corresponding  value  of  the  load 

2Nx  =  GjC  (70) 

is  referred  to  as  the  cere  shear-failure  load  in  the  classical  theory. 
When  the  bending  stiffness  of  the  faces  is  included  (Reference  8),  the 

55 


value  of  the  load  parameter  does  not  tend  to  the  asymptotic  value  but 
increases  rapidly  with  increasing  values  of  the  wavelength  parameter. 
This  case  is  described  by  Equation  (k7)  of  the  present  analysis  when 
cp  — » oo  .  The  solution  based  on  the  theory  of  Reference  5>  which  employs 
the  same  degree  of  freedom  as  the  present  development,  does  not  reveal 
the  existence  of  the  antisymmetric  wrinkling  mode  for  reasons  explained 
earlier.  That  is,  an  arbitrary,  and  obviously  eminent,  assumption  that 
the  effect  of  normal  strain  is  of  negligible  order  in  ant isymme trie 
deformation  of  the  faces  precludes  the  occurrence  of  a  possible  buck¬ 
ling  mode. 

Closed-form  expressions  for  the  critical  wavelength  parameters,  the 
associated  load  parameters,  and  the  stresses  are  shown  in  Table  II 
for  each  of  the  three  instability  modes  discussed.  It  can  be  con¬ 
cluded  from  these  results  that  the  symmetric  wrinkling  mode,  for 
practical  parameter  values,  corresponds  to  a  lower  value  of  the  load 
parameter  than  the  antisymmetric  wrinkling  mode.  If  the  converse 
were  to  be  true,  then 

2(1  +  e )  vScfT  -  0pcp  <  2  ,/qT 


or 


epep  _ 

- — —  >  /cp  (71) 

2(1  +  €)  V3  -  2 


However,  for  the  left-hand  side  of  the  inequality,  the  assumption 


P2  =  (1  +  e)  ;3qT  -  0pcp  >  1 


This  same  conclu&ior  is  evidenced  in  the  behavior  of  axially  com 
pressed,  simply  supported  columns  and  plates. 
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AXISYMMETRIC  INSTABILITY  MODES 


is  made;  that  is, 


_  1  +  0pcp 

<r  >  - z=  (72) 

(1  +  0  3 

2 

The  quantity  0pcp  =  ( 12/rr  )(f/g^)  reflects  values  from  5  to  10  for  prac¬ 
tical  honeycomb-core  materials.  If,  for  example,  e  =  0.05  ,  then  the 
conditions  Equations  (71)  and  (72)  require  that 

1  +  0pcp  <  1.82  cp  <1.11  0pcp  (73) 

and  the  interval  [1  +  0pcp  ,  1.1  0ptp]  is  nonzero  only  if  0pcp  >  9  ; 
this  appears  to  be  highly  unusual.  Thus,  the  critical  instability 
mode  is  either  the  classical  mode  or  the  symmetric  wrinkling  mode, 
according  to  the  present  theory.  It  is  also  observed  that  both 
wrinkling  modes  show  wavelengths  that  are  smaller  than  the  thick¬ 
ness  of  the  shell  and  possibly  of  the  order  of  magnitude  of  the  cell 
size  of  the  core  material.  The  extreme  situation  requires  a  treat¬ 
ment  of  the  core  as  an  aggregate  consisting  of  discrete  elements 
rather  than  an  idealized  continuum. 

As  long  as  the  theory  is  applicable,  any  intent  to  design  a  structure 
with  the  same  safety  margin  with  respect  to  different  instability 
modes  requires  that 


2(1  +  e)  3£/e2  -  0p£/e2  =  2  cp  (7*0 

or,  approximately  (c  «  1)  , 

2  2 

0prp  =  -  2  y  cp5/2  +  2  Cp  (75) 

For  each  assumed  value  of  the  curvature  parameter  £  ,  the  relation 

between  0pcp  ,  which  is  proportional  to  F/Gj  (a  characteristic 
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quantity  for  a  core  material)  and  the  core  stiffness  parameter  cp  can 
be  determined.  Figure  6  shows  a  stability  boundary  for  the  two  modes 
of  buckling.  The  ranges  of  practical  values  of  0pcp  are  indicated 
for  two  typical  core  materials. 

NONAXISYMMETRIC  BUCKLING  OF  THIN.  CIRCULAR  CYLINDRICAL  SANDWICH  SHELLS 

The  solution  to  Equation  (29)  for  the  case  of  nonaxisymmetric  buckling 
leads  to  the  buckling  criterion  given  by  Equation  (52).  The  critical 
values  of  the  wavelength  parameter  (3  and  the  corresponding  load 
parameter  k  ,  obtained  from  this  equation,  are  shown  in  Table  III. 

A  comparison  with  the  corresponding  results  for  the  case  of  axisym- 
metric  buckling  discloses  that  for  both  wrinkling  modes,  the  critical 
values  of  the  wavelength  parameter  are  l/  1  +  p  times  the  values  for 
the  axisymmetric  case.  In  the  isotropic  case  (p  =  1)  ,  the  critical 

load  parameter  is  (l  +  p)  times  as  large  for  both  wrinkling  modes. 

The  longwave  buckling  mode  (classical  mode)  also  reveals  a  higher  buck¬ 
ling  load  for  the  nonaxisymmetric  case,  as  the  positive  factor 

6  p  2e2  (1  +  p/p2) 

f  -  + - — 

p  l-v  (1  +  p) 


0  p  p  2e2 

1+p  _P(l+pp)+ 

P  l-v 

is  smaller  than  unity  if 

0p32  > 

This  condition  is  usually  satisfied,  as  all  quantities  are  nonzero  and 
positive  while  l/p2  <  1.5  ;  that  is,  the  right-hand  side  of  the  in¬ 

equality  is  negative  in  many  practical  cases. 

However,  it  should  be  noticed  that  for  shells  of  finite  length,  this 
conclusion  is  too  sweeping.  In  the  foregoing  analysis,  the  wavelength 


2^  (2  +  p  -  1/p2) 

1-V  (1  +  p)2 


(76) 
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Figure  6.  Stability  Boundary  for  Longwave  and  Symmetric 
Wrinkling  Instability  of  a  Honeycomb  Sandwich 
Cylinder  (c/R  =0.05  ,  t/c  =  0.05)  . 
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parameter  3  is  treated  as  a  continuous  variable.  When  the  cylinder 

is  assumed  to  accommodate  only  an  integer  number  m  of  halfwaves,  then 

3  =  (c/L)m  just  takes  discrete  values.  Either  one  of  the  two  discrete 

values  enclosing  the  theoretical  critical  value  of  a  continuous  variable 

3  may  correspond  to  the  lowest  value  of  the  load  parameter  k  .  This 

value  of  the  load  parameter,  of  course,  will  be  higher  than  the  minimum 

value  reached  in  case  of  a  continuous  3  .  The  magnitude  of  these 

deviations  from  the  critical  load  parameters  as  given  by  Equations  (46) 

and  (63)  can  be  ascertained  only  by  numerical  evaluation  of  a  particular 

2 

case.  However,  it  can  be  observed  that  this  effect  is  small  as  6k/6p 
is  near  zero  in  the  vicinity  of  the  critical  wavelength.  Only  in  cases 
where  the  last  term  on  the  right-hand  side  of  Equations  (46)  and  (63)  is 
small  with  respect  to  the  first  term  (0p£/e^  «  2(l+e)  )  may 

the  nonaxisymmetric  buckling  mode  yield  a  slightly  different  and  pos¬ 
sibly  lower  buckling  load.  However,  in  view  of  other  simplifications 
made,  the  axisymmetric  buckling  load  can  be  used  to  predict  the  clas¬ 
sical  buckling  mode  in  most  practical  cases. 

WEAK  CORE  SHELLS 

The  approximate  closed-form  expressions  for  the  critical  loads  and  wave¬ 
lengths  of  the  two  antisymmetric  buckling  modes  are  derived  on  the  basis 
of  assumptions  with  respect  to  the  order  of  magnitude  of  the  wavelength 
parameter  3  •  For  the  classical  mode,  3  is  assumed  to  be  much  smal¬ 

ler  than  unity.  In  order  to  assess  the  effect  of  the  core  stiffness  on 
the  critical  wavelength,  the  axisymmetric  case  is  considered  (Table  II). 
It  is  seen  that  for  the  classical  mode, 

3 

53  c  6X  3^ 

-  =  -  —  -  =  _  >  0  (77) 

6(ep)  x  5(ep)  2c* 

The  quantity  Qp  is  proportional  to  l/Gj  .  For  decreasing  core 
stiffness,  d(ep)  >  0  and,  hence,  dX  <  0  .  If  0p  approaches  the 

value  e  ( 1+e )  \JJT  >  then  the  wavelength  parameter  tends  to  infinity 
(the  wavelength  tends  to  zero)  and  the  approximation  is  no  longer  valid. 
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This  degenerate  case  is  the  classical  “shear-failure"  phenomenon,  with 
a  critical  load  parameter 


(kcr)  =  3(1  +  O2/0P  (78) 

cl 

For  the  shortwave  antisymmetric  wrinkling  mode,  the  wavelength  parameter 
P  is  assumed  to  be  of  the  order  of  magnitude  of  unity  or  larger.  To 
evaluate  the  effect  of  core  shear  stiffness  in  this  case,  it  must  be 
noticed  that  the  quantity  0pcp  is  proportional  to  F/G^  and  that  this 
ratio  does  not  vary  appreciably  for  a  given  type  of  core  material.  Then, 

8p  c  c!X  '  5"1  1  +  e 

-  =  ^  -  >  0  (79) 

Sep  X  Sep  4  '  cp  p 

But  as  dep  <  0  for  a  decreasing  core  stiffness,  dX  >  0  .  If  with  de¬ 
creasing  core  stiffness  (1  +  e)  3ep  -  0pcp  -» 0  ,  then  (p  )  -*  0 

cr  2 

(the  wavelength  increases  beyond  bounds)  and  the  assumption 
P  =  0(1)  is  violated. 

Clearly,  the  wavelengths  of  the  two  antisymmetric  buckling  modes  ap¬ 
proach  each  other  as  the  core  stiffness  decreases,  and  they  eventually 
coalesce  (Figure  5).  Which  of  the  two  critical  load  formulas  (Equa¬ 
tions  (ho)  and  (h9))  best  approximates  the  minimum  value  depends  on  the 
magnitude  of  p  for  which  the  coalescence  occurs.  The  "shear-failure" 
formula  (Equation  (78))  constitutes  a  lower  bound  in  any  case. 

The  criteria  developed  herein  are  utilized  to  predict  the  buckling  loads 
of  the  sandwich  cylinders  constructed  and  tested  as  reported  in  References 
10,  18,  and  19.  In  the  first  and  last,  axial  compression  tests  were  car¬ 
ried  out;  in  the  second,  the  cylinders  were  tested  in  bending.  With  re¬ 
gard  to  the  bending  tests,  the  cylinders  reflect  thickness -to -radius 
ratios  and  radius-to-length  ratios  sufficiently  large  to  permit  the 
use  of  the  axial  compression  buckling  criteria. 
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The  results  of  the  comparison  of  theoretical  prediction  with  experi¬ 
ment  is  presented  in  Table  IV.  The  theoretical  buckling  stresses  are 
obtained,  as  appropriate,  from  either  Equations  (1+3),  (46),  and  (1+9) 
or  their  reduced-stiffness  counterparts  developed  in  Appendix  IV  as 
Equations  (125),  (127),  and  (129).  In  view  of  -he  fact  that  all  of 
the  cylinders  tested  demonstrate  practical  construction  and  that  the 
bending  tests  were  somewhat  affected  by  end  conditions,  it  would  ap¬ 
pear  that  good  agreement  has  been  realized.  The  implication,  based 
on  the  limited  test  data,  is  that  buckling  stresses  for  sandwich  cylin¬ 
ders  with  geometric  and  material  parameters  in  the  practical  range  of 
interest  are  given  in  good  approximation  by  linear  theory  with  no  ap¬ 
parent  need  for  the  excessive  "knock-down"  factors  appearing  in  current 
design  criteria.  The  ratio  of  experimental  and  theoretically  predicted 
critical  stress,  t]  ,  based  on  the  data  of  References  10,  18,  and  19 
and  the  present  theory  (see  Table  IV),  is  shown  in  Figure  7.  For  com¬ 
parison,  design  recommendations  applicable  to  general  instability  of 
initially  imperfect,  linear  elastic  shells  (see  References  20  and  21) 
are  included  in  Figure  7- 
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Figure  7.  Ratio  of  Experimental  and  Theoretical  Buckling  Stresses  for  Perfect 
Shells  (Table  IV)  Compared  with  Design  Recommendations  for  General 
Instability  of  Linear  Elastic,  Imperfect  Shells  (References  20  and  21). 


CONCLUS IONS 


This  report  presents  the  governing  equations  and  boundary  conditions 
for  moderately  large  deflections  of  a  sandwich  shell  of  arbitrary 
shape.  Equilibrium  equations  and  boundary  conditions  for  both  the 
thin-face  layers  and  the  thick-core  layers  are  consistently  derived 
from  the  condition  for  a  stationary  value  of  the  total  potential. 

For  the  core,  the  face  parallel  stresses  are  neglected;  but  the 
transverse  strain  is  included,  and,  unlike  previous  theories  em¬ 
ploying  the  same  degrees  of  freedom,  the  relevant  strain  displace¬ 
ment  relations  for  the  core  are  linearized.  The  stiffness  parameters 
of  the  core  are  assumed  to  vary  linearly  through  the  thickness,  as  is 
the  case  in  a  curved-core  layer  of  cell  structure.  In  Reference  5> 
not  only  linear  terms  in  the  normal  coordinate  are  retained  but  also 
higher  order  terms,  while  the  stiffness  parameters  are  assumed  to  be 
constant.  The  present  approach  is  believed  to  provide  a  more  consis¬ 
tent  treatment  of  the  core.  These  assumptions  allow  the  three  equi¬ 
librium  equations  of  the  core  equation  to  be  integrated  with  respect 
to  the  normal  coordinate;  all  loads  transferred  by  the  core  can  be 
expressed  in  terms  of  the  components  of  the  face  displacements  and 
two  core-shear  angles,  while  two  compatibility  equations  are  obtained 
as  auxiliary  conditions. 

As  a  final  set,  three  equilibrium  equations  and  four  boundary  condi¬ 
tions  for  the  edge  of  each  of  the  face  layers  are  presented  in  terms 
of  the  face-displacement  components.  They  are  of  the  type  usually 
obtained  on  the  basis  of  the  Kirchhoff-Love  assumptions  for  a  thin 
shell. 

In  addition,  two  compatibility  conditions  and  one  boundary  condition 
for  the  core  are  given  also  in  terms  of  the  face  displacements.  The 
eight  governing  equations  and  the  nine  boundary  conditions,  then,  are 
functions  of  two  surface  coordinates  only. 
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Linearization  leads  to  a  set  of  equations  with  a  symmetric  operator 
matrix  as  is  appropriate  for  a  linear  elastic  system.  The  order  of 
a  single  governing  equation,  as  derived  from  the  operator  matrix  of 
the  full  set,  is  lQ,  with  9  boundary  conditions  which  can  be  speci¬ 
fied  at  an  edge  in  terms  of  either  generalized  force  or  displacement. 

The  specification  of  a  consistent  set  of  boundary  conditions  is  nec¬ 
essary  not  only  for  the  derivation  of  an  exact  or  approximate  solution 
to  a  problem  but  also  for  a  meaningful  comparison  of  theoretical  pre¬ 
dictions  and  the  results  of  appropriate  tests.  Although  the  nine  sets 
of  boundary  conditions  occur  in  a  discouragingly  large  number  of  pos¬ 
sible  combinations,  practical  designs  will  limit  the  number  of  com¬ 
binations  and  should,  in  addition,  stimulate  the  analysis  of  important 
boundary  situations  other  than  the  very  limited  set  that  is  usually 
considered. 

The  unified  theory  is  used  to  estimate  the  range  of  applicability  of 
classical  theories  or  simplified  versions  of  the  present  theory.  In 
an  example,  the  critical  loads  for  an  axially  compressed  circular 
cylinder  are  determined. 

In  general,  three  instability  modes  are  possible: 

1.  An  instability  pattern  with  a  wavelength  of  several  times  the 
shell  thickness  and  negligible  transverse  strain.  This  mode 
corresponds  to  the  results  obtained  from  the  classical  sand¬ 
wich  theory,  ignoring  core  normal  strain. 

2.  A  face-wrinkling  pattern  of  very  short  wavelength  (of  the 
order  of  the  shell  thickness)  and  characterized  by  essen¬ 
tially  similar  face  displacements  (antisymmetric  wrinkling). 

3.  A  face -wrinkling  mode  of  comparable  wavelength  but  face  dis¬ 
placements  that  are  essentially  symmetric  with  respect  to  the 
middle  surface  of  the  shell  (symmetric  wrinkling). 

By  means  of  a  parameter  study,  approximate  expressions  for  the  critical 
loads  and  wavelengths  of  all  three  instability  modes  are  derived.  The 
limits  of  applicability  are  indicated.  The  study  also  reveals  that  the 
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symmetric  wrinkling  mode  will,  for  most  practical  sandwiches,  lead  to 
the  lowest  wrinkling  load.  In  Reference  3,  the  opposite  conclusion  is 
reached;  however,  it  appears  that  the  results  of  numerical  analysis 
are  not  consistent  with  the  basic  theory. 

The  parameter  analysis  further  discloses  that,  with  decreasing  core 
stiffness,  the  wavelengths  of  the  two  antisymmetric  modes  (wrinkling 
and  overall  instability)  approach  each  other  and  finally  coincide. 

The  coalescence  of  critical  wavelengths  can  occur  near  the  overall 
instability  wavelengths  as  well  as  near  the  wrinkling  wavelength. 

Only  in  the  first  case  is  the  classical  "shear  failure"  formula  a 
good  approximation. 

The  critical  loads  derived  from  a  nonaxlsyrnnetric  deformation  pattern 
may  be  lower  in  case  of  the  longwave  deformation,  but  it  is  shown 
that  the  difference  will  be  small. 

A  comparison  with  limited  experiments  suggests  that  critical  loads 
derived  from  a  linear  analysis  of  axially  compressed  sandwich  cylin¬ 
ders  should  be  used  to  predict  the  buckling  loads  of  such  cylinders. 

The  present  analysis  shows  also  that  a  consideration  of  the  discrete 
core  structure  may  be  necessary  to  account  for  existing  discrepancies 
between  linear  theory  and  experiment.  However,  unlike  the  case  of 
homogeneous,  isotropic  shells,  the  discrepancies  that  do  exist  appear 
to  be  more  a  function  of  core  deformation  prior  to  buckling  and  in¬ 
elastic  behavior  at  buckling  rather  than  the  result  of  overall  im¬ 
perfections  in  the  sandwich  shell  wall. 
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APPENDIX  I 

VARIATIONAL  DERIVATION  OF  THE  GOVERNING  EQUATIONS  FOR  FACE  SHEETS 


The  face  sheets  are  treated  as  thin,  shallow  shells  of  linearly  elastic 
isotropic  material.  The  governing  equations  are  derived  from  the  con¬ 
dition  for  a  stationary  value  of  the  total  potential.  The  total  poten¬ 
tial  is  defined  as  the  sum  of  the  strain  energy  function  and  a  potential 
function  of  the  applied  loads,  as  appropriate. 


THE  STRAIN  ENERGY  FUNCTION 

The  appropriate  strain  energy  density  expression  for  linearly  elastic 

7 

thin  shells,  consistent  with  the  Kirchhoff-Love  assumptions,  is: 


+  2V£1£2 
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For  the  case  of  small  strains  and  moderate  out-of-plane  rotations,  the 
inplane  rotation  about  the  normal,  determined  by  the  quantity 
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is  neglected  with  respect  to  the  rotations  of  the  normal  to  the  middle 
surface  given  by 
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The  middle  surface  strains  are  then  approximated  by 


The  geometry  is  shown  in  Figure  1;  the  coordinate  directions  1  and  2 
are  directions  of  the  principal  curvature;  the  principal  curvatures 
are  l/pj  and  l/pg  ,  respectively.  The  changes  of  curvature  and 
twist  of  the  deformed  middle  surface  are 


The  total  strain  energy  of  a  face  sheet  is  then 

Uf  =  JJ [w^c^  ]  d^d^  (i  =  1  or  2)  (84) 

*  D  i 
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POTENTIAL  OF  APPLIED  LOADS 


The  boundary  loads  on  a  face  sheet  are  indicated  in  Figure  2.  Upper 
case  letters  denote  loads  per  unit  length;  lowercase  letters  denote 
loads  per  unit  area. 

For  the  variational  problem,  as  the  equivalent  of  the  principle  of 
virtual  displacements  for  an  equilibrium  system,  a  potential  of  ap¬ 
plied  loads  V  is  appropriately  defined  as: 


V, 
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in  which 
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The  total  potential  function  is  then  defined  by 


Ti  =  Ui  +  Vi 


(86) 


It  should  be  noted  that  the  neglect  of  the  inplane  rotations  and  the 
product  terms  in  the  potential  expression  involving  the  out-of-plane 
rotations  leads  to  Euler  equations  of  a  first  variation  that  are  incor¬ 
rect  in  terms  involving  rotation  components.  Adjustments  will  be  made 
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later  to  obtain  the  simpler  form  of  equilibrium  equations  for  coordinate 
directions  on  the  deformed  surface  rather  than  for  the  original  coor¬ 
dinate  directions. 


Variation  of  the  Total  Potential 

For  an  equilibrium  system,  the  first  variation  of  the  total  potential 
vanishes  for  any  set  of  variations  of  the  argument  functions  compatible 
with  imposed  boundary  conditions.  Carrying  out  the  first  variation  and 
integrating  by  parts  integrals  containing  derivatives  of  the  argument 
functions  u^  ,  v^  ,  leads  to 
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The  potential  term  involving  the  twisting  moment  has  been  integrated 

by  parts  once  with  respect  to  the  circumferential  coordinates  S  .  The 
contour  is  assumed  to  be  smooth;  if  the  contour  contains  singularities 
(corners),  the  integration  by  parts  also  gives  rise  to  terms  that  are 
the  product  of  a  variation  in  the  lateral  displacement  w  and  a  con¬ 
centrated  generalized  force  ("corner  reaction").  It  is  well  known  that 
this  identity  implies  that  each  of  the  integrals  must  vanish  identically. 
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As  the  variations  Su  ,  Sv  ,  and  5w  are  completely  arbitrary  within 
the  domain  of  integration  S  ,  the  condition  ST  =  0  leads  to  the 
following  set  of  Euler  equations  and  boundary  conditions.  Appropriate 
quantities  have  to  be  assigned  a  subscript  1  or  2  .  Where  a  double 
sign  appears,  the  top  sign  is  to  be  used  in  connection  with  the  sub¬ 
script  1  ;  the  lower,  with  the  subscript  2  .  In  S  ,  then, 
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c(€i  +  V62)  cuc21  +  c  —  r12cu  +  c(€2  +  V€j)  c21cu  +  c  —7 12c21 


c-  c  c 

-  D(l-v)  K12  —  -  D(<1+VK2)  11  —  -  D( Kg+VKj) 

Pi 


c21cll 


Coi  M  M 

D(l-V)  K  -  +  N  -  -S  +  — 

p2  T  Ps 


=  0  or  5dsl  =  0  (89a) 
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l-v 

C(cl  +  '"-2)  c12c21  +  c  “  >'l2c12cll  +  C(c2  *  vcl>  C22CU 


1”V  C,,C 


+  °  2  712C22C21  "  D(1‘v)  *12  p 


2 

iri2 


1 


-  D(<2  +  VKX)  C22c 1 1  "  D(1_v)  K 


°21C22 


12 


d(kx  +  \*2) 

p, 


M  M 

+  N  -  -2  +  —  =0  or  6d  .  =  0 

n  „  T  ni 

Pn  T 


C(c 


(89b) 


(1  u  \  l-v  /  1  V  V 


+  C(g2  +  ve 


/I  vv  1-V/l  u  V 

j)  —  w,.  -  —  c  +c  —  rig!  —  w,ft - Jcn 

1  \a  h  P2'  2  'ai  1  Pi* 


a 


d„  |d(k,  +  -/K0)  c„  +  D(l-v)  itjjC..! 


'll  |  "'"1  T  2’  21 


11 1 


*  d12  Jd(ks  +  vKj)  cn  +  D(l-v)  Kj2c2l  j 


’8. 


ala2 


J  a2  d(ki  +  *2)  [  ,  -  d(k2  +  vkx)  a2>ii 


S  ■><>-v>«18j,|  +  b:»-v)  *isPx,i 


21 


aia2 


jaj  d(k2  +  VKX)  J  -  d(kx  +  VKg) 


(Cont*d) 
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+  ja2D(l-v)K12|  +d(1-v)k12o.2(6i 


11 


t  t  1 

+  :  v?i +:  vn  -N*  +~Vi 

2  2  a  *  8 

s 


=  0 


or  =  0 


(89c) 


Jd(Kj  +  VKg)  c21  +  D(l-v)  j  d21  +  Jd (k2  +  vKj)  cn 


+  D(l-v)  <12cgl J  d22  4  Mn 


=  0 


or  8w  .  =0 

l'sn 


(89d) 


The  three  Euler  equations  for  the  variations  5u  ,  5v  ,  and  8w 
(88a, b,c)  are  recognized  as  the  equilibrium  equations  for  components 
of  forces  In  the  three  coordinate  directions.  The  last  Euler  equation 
has  been  rewritten,  using  the  first  two. 


The  four  sets  of  boundary  conditions  (89a,b,c,d)  are  obtained  from  the 
boundary  Integrals  by  use  of  suitable  transformations  from  the 
coordinate  system  into  the  £s;&n  system.  The  transformation  symbols 
c  __  and  d  are  defined  by 


c 


li 


d 


11 


Hi 

*1 

ai 

a 

as  ’ 

I  « 

1  ro 

s 

8 

n 

n 

a 

dlg 

a 

8 

d  =  — 

s 

Q1 

^1 

12  a2 

a52 

etc. 


etc. 


(90) 
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APPENDIX  II 

THE  GOVERNING  EQUATIONS  FOR  A  CORE  MODEL 


The  analysis  of  the  ''ore  Is  based  on  the  assumption  that  the  inplane 
shear  stiffness  ano  <?*tensional  stiffness  of  the  core  are  negligible 
with  respect  to  the  transverse  shear  and  extensional  stiffness  quanti¬ 
ties. 

Further,  if  the  rotation  about  the  normal  to  the  1-2  coordinate  surface 
is  considered  to  be  much  smaller  than  the  rotation  about  tangents  to 
the  same  surface  (consistent  with  the  assumption  made  for  the  face 
sheets),  then  the  essential  strain-displacement  relations  reduce  to 

7 13  = 

723  =  -^a2,C  +  ^w'i2 

12  12 
3  “  W,C  2  U’i  2 

For  a  core  with  a  cell  structure  consisting  of  walls  perpendicular  to 
the  t;  =  0  coordinate  surface,  the  last  two  terms  in  the  expression 
for  the  transverse  strain  correspond  to  rotation  of  cell  walls  out  of 
their  planes.  If  the  squares  of  these  rotations  are  not  negligible  with 
respect  to  the  linear  term  w,^  ,  then  the  effects  of  the  stresses 

and  strains  in  the  cell  walls,  that  correspond  to  a  state  of  plane 
stress  and  are  neglected  a  priori,  would  certainly  have  to  be  retained 
in  the  equilibrium  equations. 

Consistent  with  the  initial  neglect,  therefore,  the  strain-displacement 
relations  of  the  core  are  linearized;  that  is,  the  underlined  terms  in 
Equation  (91°)  are  omitted. 


(91a) 

(91b) 

(91c) 
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For  a  core  of  cell  structure,  the  appropriate  stiffness  quantities  are 


F^^Ofe  =  F(  s 

Gl^)a2  =  Gl^°)“2  S  Gl“2 


G2(UQ!1  =  GgCO)^  s  GgO^ 


(92) 


The  strain  energy  of  the  core  is 


U  = 
c 


f/f[s  ♦iVO  +  I  g2(U  rj3 


a^d^dlgd  i  (93) 


and  the  potential  of  the  applied  loads  is 

c/2 

I  |(qxu  +  q2v  +  q?w)  O^Otg  | 

s 

The  total  potential  of  the  core  is 


Vc  =  -  *  V  +  q5w)  alaS  1  d5ld62  • 

s  L  ”c'2  as  -c/2 


q  «*gd£sd£ 


(9b) 


T  =  U  +  V„ 
c  c  c 


(95) 


For  a  stationary  character  of  the  total  potential,  the  first  variation 
must  vanish  Identically;  thus, 


5T  = 


Iff 


C  JJJ  5u  r  ^alGla2713^  “  Cla2713°l,c 
R  ^ 


fff  8v  ['  ^a2G2al723),t  "  G2ai723a2,C 

R  = 


d^dftgdg 


dexdS2dC 


fff  Sw  -  BttjOg  w,cc  -  (Gxa2715)  -  (GgO ^23) 


d^dtgdt 
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d*ld*2 


ft 

ft 


c/2 

+  *p  r  5w 

8S  -c/2 


(  a2  1 

,c/2 

Gl-^3-"l 

2  '  ■ 

CVJ 

o 

1 

11 

*_n 

(  “l  ) 

c/2 

jG2~723 

5-c/2 

(  5152  ) ' 

,c/2 

r  ’Vb 

(  axa2  b  ^JJ 

CVJ 

o 

i 

ii 

„  ^13«2«2  +  G2 
a2 

ai 

“  7  23 
al 

d|,d| 


1  *2 


d*ld*2 


d£ 

(96) 


As  the  variations  6u  ,  8v  ,  and  5w  are  arbitrary  everywhere  Inside 
the  domain  of  Integration  and  In  the  interfaces,  the  following  Euler 
equations  and  boundary  conditions  result: 


(alGlVl3>  +  GlVl3al,t  ’  ~~  (alr15)„  =  ° 


(a2G2aiy23\t.  +  G2ai723a2,C 


+ 


a. 


G„a. 

— '  ^a2723^ 
a2  2  25  't 


-  (a2Gl713^  +  ^alG2723^  1 

WjOfe  2  1  w  'lx  2  25  'tj 


(in  R) 


«  0 


=  0 


(97) 


a. 


Gi  ”  713  =  q 


a. 


a. 


'23  -  "2 


(in  S  ,  £  =  ±  c/2) 


(98) 
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a.a 

F  w, 

ala2 


■ 


G1  2  r13a2d52  +  G2  1  r23ald5l 

a2  1 


q  a  d  or  6w  =  0 

^  i  i 


(on  dS  ,  |U  <  c/2)  (99) 


The  Euler  equation*  can  be  Integrated.  The  first  two  equations  yield 


713(?)o^(C)  =  713(0)0®(0)  = 
725(C)<^(C)  -  723(0)o|(0)  =  7  2^ 


Further,  as  0^(1;)  *  a^(0)(l  +  £/pt)  ,  these  results  can  be  written  aa 

r15(t)  -  rt  (1  -  2C/ox)  +  o«s/p=) 

r2J(?)  -  ra  (l  -  «/p2)  *  o{?2/P|)  (loo) 

The  shear-stress  distribution  across  the  thickness  of  the  core  model  la 
completely  determined  by  the  equilibrium  equations. 


The  last  Euler  equation  can  be  Integrated  twice  with  respect  to  £  . 

With  the  use  of  Equation  (100),  there  results 


/?  t,2  \  c  _  ( f  c2  \  c2  _ 

"(SwU  )  =  [  —  -  —  )  “Z-Z-  (a27i)  +  (  —  •  —  )  izr  (<y2) 

1  2  V 3Pj  2  /  FV2  '51  '5P2  2  '  ^2°!  'h 


+  C  +h(l1,l2) 


07 


o  o 

Terns  of  Che  order  £  /p^  have  conaisCenCly  been  neglected  with  reapect 
to  unity. 

The  two  functions  £( )  and  h(|j,S2)  can  be  determined  by  making 
use  of  the  "boundary  conditions" 


wfSjiSgjc/S)  =  Wj 

w(5j  ,  i2  ,  -c/2)  =  w2 


(101) 


These  relations  actually  express  continuity  of  the  lateral  displacement 
w  of  the  interfaces,  as  the  lateral  strain  in  the  faces  la  neglected. 
By  use  of  Equation  (101),  there  results 


+  w_ 


2 


8  12Pj 


8  12Pg 


and,  consequently. 


-  w„ 


i!+ii 

2  3p, 


i!+ii 

2  Jp, 


~t*t  (Vj) 
rai°2  ,{1 


—  — 


t. 


Px  12Pj  /  F ttjCfg 


(Vi) 


,2  _2 


pg  12p2  /  FttgSj 


(®1^2^ 

L  £  n 


(102) 


(103) 
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From  the  strain  displacement  relations  (Equations  (91a)  and  (91b)),  it 
is  seen  that 


1' 


't 


1  1 

y13  "~2  w*  | 

ax  D  cq  5l 


a.  r®  '  o?.V'h 


Integration  of  tiiese  relations  with  respect  to  £  ,  making  use  of  pre¬ 

vious  results  and  the  continuity  conditions, 


-f) 

vfSj.tj,  -|) 


2 

v 


) 

) 


results  In  two  compatibility  conditions  in  terms  of  the  components  of 
the  face  displacements  and  the  core  shear  functions  7^  and  7^  ; 

namely. 


0  (105) 


89 


V  /i  )-  v2(i  +  —  )  +  r-»! 

\  2p J  '  2Pg/  a2  W 


c  +  t1  C  + 

2  6Po 


1  c  +  t  c  +  3t  c 

+  —  w  -  + - 

a2  ,52  2  6p 


-  7gc 


12Ftt, 


G1  -  Go  _ 

ITT-  (<V,)  +  ~  (a-,7p) 

axa2  a2ax  d  ' 


=  0  (106) 


2  >  §r 


The  boundary  conditions  (Equation  (98))  allow  elimination  of  the  inter¬ 
face  loads  q^  ,  q^  ,  and  q^  from  the  Euler  equations  of  the  faces; 
then  by  use  of  Equations  (100)  and  (103),  resulting  from  integration 
of  the  equilibrium  equations  of  the  core,  the  only  unknown  functions 
remaining  are  the  six  face  displacement  components  and  the  two  core 
shear  functions.  Each  of  these  eight  functions  depends  on  the  surface 
coordinates  and  only. 

In  addition  to  these  six  equilibrium  equations  for  the  faces,  the  two 
compatibility  conditions  (Equations  (105)  and  (106))  have  to  be  satisfied 
to  guarantee  a  unique  solution.  The  last  boundary  condition  for  the  core, 
Equation  (99),  also  has  to  be  satisfied  in  addition  to  the  eight  boundary 
conditions  obtained  for  the  two  faceplates. 

If  it  is  assumed  that  the  boundary  load  q  is  distributed  over  the 
depth  of  the  core  in  the  same  way  as  the  internal  shear  stresses, 


q*(S)  -  q*(o)  (l  -  2?/p„)  -  q  * (l  -  2C/p„) 


then  the  ninth  boundary  condition  can  be  written  in  final  form  as 


Vlc21 


*  Vscn 


-  * 
+  q 


0  or 
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Bi  (x-;rM  1+;r) 

fl  R 


C  G1  -  Go  - 

+  —  ZTT  (Vi)  +  (<Vl) 

6f  (  c^a  ’  ^  a  a  '  t 


=  o  (107) 


According  to  Salnt-Venant's  principle,  replacement  of  the  actual  boundary 


q  (O 


by  the  statically  equivalent  load  q  ca  will  affect  the  stress  dls- 
tribution  in  a  limited  edge  zone;  it  also  affects  the  plane-stress 
components.  However,  as  these  stresses  have  been  ignored  in  the 
present  theory,  Equation  (107)  is  the  appropriate  boundary  condition 
of  the  core.  If  the  plane-stress  components  are  Important,  then  a 
more  complete  treatment  of  the  core  requires  consideration  of  other 
than  shear  loads  at  the  boundary. 
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APPENDIX  III 


BUCKLING  EQUATIONS  FOR  AN  AXIALLY  COMPRESSED  CIRCULAR 
_ CYLINDRICAL  SANDWICH  CONSTRUCTION _ 


For  an  axially  loaded  cylinder,  the  shallow-shell  equations  can  be  used 
with  great  confidence  to  determine  buckling  or  postbuckling  behavior. 
The  following  substitutions  are  made  in  order  to  specialize  Equations 
(19)  and  (20)  for  this  case  (see  Figure  3).  For  the  reference  surface, 


h  -  x  “l  "  1  P1  ■  "  -  pn 


'  <p  “2  ■  R  *  pa  -  p. 


For  face  1  (outer  face): 


0^  =  1  p^  *  00 


a.  =  R.  =  pr 


For  face  2  (inner  face): 


0^=1  Px  =  “ 


a2  R2  P2 


At  the  boundary,  x  =  0  : 


(108) 


K  =  -  <P  a«  = 


Ri  =  ps 


d  =  -  v  N  =  -  N  M  =  -  M 


s  cp 


<P 


I  =  x  a  = 


n 


1  p  =  00 

n 


d  =  u  N  =  N  M  =  M 

n  n  x  n  x 
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cn  -  0  -  css  r  =  0  Nz  =  \ 


c12  "  1  "  C21 


dH  -  0  -  d22 


d12“'1“-d21 


(109) 


At  the  boundary,  x  =  L  : 


iB  -  <P 


a  = 


Ri  =  P8 


d  =  v 
s 


N  =  N 

s  cp 


M  =  M 

8  (p 


«n  =  * 


a  = 

n 


1  pn  =  “ 


d  =  -  u 
n 


N  =  -  N  M  =  -  M 
n  x  n  x 


Cll=  0  =  C22 


T  =  0  N  =  -  N, 

z  t 


C12  ~  -1  =  ”C21 


dir  0  -  d22 


d12=  1  "  ‘  d21 


(no) 


Compatibility  conditions: 


Uj  -  u2  -  cyj  + 


c 

12F 


Gl7l,x  +  G272,^ 


c  +  t  c  +  t 

+  - w,  + -  W-,  =0 

0  1>X  O  2>x 


(111a) 
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/  c  \  / 

c  \  c5 

1  +  —  J-  C7  +  - 

2R  /  12 F 

G,7,  +  G-7-  — 

U,x  2'2,cp 

»<p 


°‘n, . 

I  s  6r  rl«» 


rc  +  t  c  +  3t  C 
+  I - -  + - — 


f 


6r 


K;  ■  0 


(mb) 


Equilibrium  equations  for  faces  1  and  2 


1-V 

Cl(,:l  +  v£S),x  +  C1  y  ?12,¥  T  0 


m  f1  t  i)  “  °  (nic) 

V£2  +  v£lU  +  Ci  T  r12,x  *  Vs  (l  T ; )  -  0  <llM> 

V'l  +  "eKx  +  Dl(1-',>  *12,,.  +  Dl(*2  +  '“I),,  +  Dl(l-V>  *12,, 

J  fX 


-  Ci(e1  +  VGg  Kj  -  Ci(€g  +  VP 


/  1  \  1_v 
l>  r2'^)'8C‘T 

*  [KK«*K) 


712«12 


fg7  _ 

F  72/<P 


G272,9  =  0 


(HI®) 


Boundary  conditions  at  x  =  0,L  : 

C1<£1  *  V£P>  ■  Nxi  = 


5u1  =  0 


9b 


0 


or 


(112a) 


-  +  v*2)  -  M  £  =  0  or  6w1 


(llSd) 


Under  the  shallow- she  11  assumptions,  the  strains  and  changes  of  the 
curvature  are  simply 
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If  it  is  assumed  that  the  cylinder  deforms  uniformly  under  a  uniformly 
applied  axial  compressive  load  up  to  a  critical  value  of  the  load  (that 
is,  no  edge  loads  are  present  other  than  the  compressive  force  Nxi  ), 

then  the  subcritical  state  of  deformation  is  characterized  by  the  uni- 

«#■ 

form  lateral  displacements  w^  and  Wg  ,  the  absence  of  shear  deforma¬ 
tions  y ^  and  y g  of  the  core,  and  the  circumferential  displacements 
v^  and  Vg  .  Also,  in  the  present  notation,  the  changes  of  curvature 
as  well  as  derivatives  in  the  ep-direction  are  zero. 


From  (111a)  then,  it  follows  that  u^  =  u^  =  u  ,  Equations  (lllb,d) 
and  the  boundary  conditions  (112b,c,d,e)  are  satisfied  identically. 
Equations  (111c)  and  the  boundary  conditions  (112a)  show  that 


* 

NXl  =  Ct(€1+V€2)  =  Ct(u*x  +  V^j 


while  Equation  (llle)  becomes 


r  * 

W1  * 
—  +  VU. 


1  L  l 


w, 


+  F 


*  * 

r 

2 


(c  \  w.  -  w, 
1  --  - — 1 

2R  /  c 


_2  * 
+  vu,. 


-  F 


(i*-) 

\  2R  /  c 


=  0 


#  * 

t 

2 


(11M 


(115) 


Subcritical  deformations  are  marked  by  an  asterisk. 


Multiplication  of  the  first  of  Equation  (115)  by  ^  and  of  the  second 
by  Rg  shows,  after  addition  of  the  two  results,  that  the  total  cir¬ 
cumferential  force  vanishes;  that  is, 
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Finally,  solution  of  Equations  (114)  and  (115)  yields 

c/R 

cT  cTc 

1  +  —  (1  -  c/2R)  +  (1  +  c/2R) 

Cx  FR 

c/R 

C  C  c 

1  +  —  (1  +  c/2R)  +  (1  -  c/2R) 

C2  Fk 

This  determines  the  ratio  of  and  necessary  to  achieve  a 

uniform  prebuckling  deformation. 

For  a  thin  shell  c/R  «  1  and,  in  this  case,  neglecting  terms  of  the 
order  c/R  gives 


Nxi  =  u*x  (H7) 

Poisson's  ratio  has  been  taken  to  be  the  same  for  both  faces.  For 
simplicity,  attention  may  be  restricted  to  the  case  of  completely 
similar  facings,  (Et^  =  (Et)g  .  Then,  simply, 


xl 


=  N 


x2 


N 


(118) 


If  the  prebuckling  deformations  are  subtracted  from  the  total  deforma¬ 
tions  and  the  buckling  deformations  so  obtained  are  substituted  into 
the  strain  and  curvature  expressions  (Equation  (115)),  Equations  (ill) 
and  (112)  describe  the  postbuckling  behavior  of  the  thin  sandwich  shell 
under  the  assumptions  made. 

For  a  determination  of  the  critical  load,  products  of  buckling  strains 
and/or  stresses  can  be  omitted  and  the  "buckling  equations"  for  a  thin 


M 


I 


i 


i 
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sandwich  cylinder  with  similar  facings  are  obtained  as: 


”  ("I,,  *  ^  }|. 


+  c 


(1-v) 


x  2 

(1-v) 


U-  +  V, 

l,cp  l,x 


>9 


-  Gl7l  =  0  (119a) 


u_  +  v  |  +  —  I  +  C  - 

2,x  \  2,cp  »  /  x  2 


U2,cp  +  V2,xJ^  +  Vl  =  0  (119b> 


w. 


(l-v) 


v,  +  —  +  vu,  +  C 

^  R  l»*  ],„ 


u,  + 
l,cp 


V*'*L 


G272  =  0  (119c) 


w„ 


V2,<f  +  f  +  W2,x 


(l-v) 

+  c - 

,cp  2 


U2,cp  + 


V2  x 

2’*J 


+  G2r2  =  0  (I19d) 


DV  w  -  - 

2 


Glrl,x  +  Va,* 


+  F 


w  -  w  c 

~  ^  (Gl7l,x  +  G272,^ 


c 

+  - 
R 


w. 


ivi  +  —  +  vu, 

[  1,<P  R  M 


+NxWl,xx=0  ^119e) 


4 

DV  w_ 


G,y,  +  G_7_ 

11,*  2  2,cp 


-  F 


wi  -  w2  c 


c  2F 


+  -  ^Glrl,x  +  W 


c 

+  - 
R 


V2,cp  +  -  +  VU2,x 


Nxw2,xx  *  0  <U9f> 


t  +  c 


U2  + 


t  +  c 


c 
12F 

,3 


(»1  +  w2)>x  -  rxc  +  —  (Vi,x  +  =  0  (119g) 


V,  - 


v2  + 


(W1  +  "2),9  '  V  +  ~  (Gl7l,x  -  *  0  (U9h> 


Appropriate  boundary  conditions  for  the  problem  under  consideration  are 
to  be  chosen  from  the  following  set  of  pairs: 
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At  x  =  0,L  : 


V  *  v<vl,„  *  Wl/Rl) 


=  0  or  Uj  =  0 


( 120a ) 


u2,x  +  v(v2,„  +  W 


=  0  or  Ug  =  0 


(120b) 


(l-v) 


u,  +  v, 
1,<P  l,x 


=  0  or  Vj  =  0 


(120c) 


U-v) 


u2,cp  +  V2,x 


0  or  Vg  =  0 


(120d) 


Wl,xxx  +  Wl,qxpx 


- 1  Vi  -  Vi,x  =  0  or  "l  ■ 0  (l80e) 


w2,xxx  +  f2'v)  w2(Wx 


-  |  GjXj  -  V'e.x  '  0  01  ”2  *  °  <120f> 


W,  +  VW, 

l,xx  l,cpcp 


=  0  or  w,  =  0 

l,x 


(1208) 


W-  +  VW- 

C,xx  2,qxp 


=  0  or  v0  =  0 


(I20h) 


CG 


l'l  '  0  or  wi  +  w2  +  fr  <Vi,x  *  =  0  (1201) 
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APPENDIX  IV 


SIONS  FOR  C 


The  lose  of  stability  of  a  sandwich  cylinder  under  uniform  axial  com¬ 
pression  may  occur  at  stress  levels  beyond  the  elastic  limit.  In  this 
case,  the  change  in  atrain  energy  due  to  a  small  perturbation  of  the 
prestressed  state  has  to  be  determined,  taking  into  account  the  sign 
of  the  changes  in  strain.  Stress  increments  and  decrements  are  re¬ 
lated  to  changes  in  strains  by  different  laws.  This  results  in  modi¬ 
fied  expressions  for  the  bending  stiffness  and  extensional  stiffness 
of  the  faceplates  and  presumably  also  for  the  core  stiffness. 


For  the  behavior  of  the  faceplates  under  e  prestress  beyond  the  elastic 
limit,  the  following  generalization  of  Hooke's  law  is  used: 


#  # 

where  E  =  E  or  Efc  and  |i  =  v  or  n  according  as  a  stress  inten¬ 
sity  decreases  or  increases.  E  and  Et  are  Young's  modulus  and  the 
tangent  modulus  as  obtained  from  a  uniaxial  stress-strain  curve  (pos¬ 
sibly  in  a  tensile  test),  and  v  and  n  are  Poisson's  ratios  for  the 
elastic  and  the  Inelastic  range;  plastic  deformations  are  assumed  to  be 
volume  preserving  (4  =  l/2)  .  The  stress  intensity  <r  is  defined  here 
as  the  square  root  of  the  second  invariant  of  the  stress  deviator  ten¬ 
sor  J2  .  In  the  case  of  the  usual  plane  stress  approximation  for 
face  sheets. 


—  (o* 

Y?  x 


o-  cr 
x  cp 


+  )‘/2 

xcp' 
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i 

1 


so 


8a  =  — 
6a 


(a  -  a  )  (8a  -  8a  )  +  6t  8t 

'  x  qr  '  x  qr  xcp  xcp 


For  a  perfect  cylinder  in  the  prebuckled  state,  a^  =  -=  0 


8a  - 


a  (8a  -  8a  ) 

X  X  CD 

2  Kl  vT 


Hence, 

(122) 


SYMMETRIC  WRINKLING  MODE 

For  this  mode,  =  -  Wg  ,  while  the  other  displacement  components  and 
the  core  shear  angle  are  essentially  zero.  The  only  perturbations  of 
the  prestrained  state  are  a  bending  strain  component  8e  =  -  zw,  and 

X  XX 

a  circumferential  strain  component  8e^  =  w/R  .  Then, 


8a  - - - -  (8e  -  5e  ) 

2  (1  +  /)  VT  X  * 


But 


Kl  /  1%!  ■ 


2z 

2 

Rtrr 

2z 

t 

2X& 

t 

e2/yT  »  i 


except  near  the  middle  surface  (z  =  0)  ,  so  the  sign  of  8a  is  the 
opposite  of  the  sign  of  8cx  and  the  face  bending  stiffness  is  taken 
as 


*  -  +  t. 

D  ~  E  “ 
12 


(125) 


—  * 


where  the  reduced  modulus  E  is  the  well  known  double  modulus: 

i  rs— 1\-2 


-  *  . 
E  =  4 


^l-v2 


Ue  e. 


w 


(ve  +  v\7 
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According  to  Equation  (40),  the  wrinkling  load  is  determined  by 


k 


<T  tC 
X 

2  * 
7T  D 


2  *,  2 
(3  +  cp  /r 


=  P2  + 


*  3 
2F 

k*p 
7T  D  p 


(124) 


In  this  expression,  F  is  a  reduced  core  modulus,  As  the  core  is 
assumed  to  be  stress-free  up  to  the  point  of  instability,  F  should 
be  equal  to  the  elastic  value  F  .  Actually,  the  core  is  attached 
to  the  faces,  and  the  outer  core  layers  will  undergo  the  same  inplane 
strains  as  face  layers  that  are  strained  beyond  their  elastic  limit. 
In  addition,  buckling  of  cell  walls  may  be  a  source  of  deterioration 
of  the  core  stiffness. 


An  upper  limit  for  the  wrinkling  stress  is  found  by  taking  F  =  F  ; 

—  _ 

a  low  value  results  if  it  is  assumed  that  F  =  F»E  /E  .  In  both 

•X-  *  ** 

cases,  both  D  and  F  are  completely  determined  by  the  face  stress 
cr^  and  will  be  stationary  if  crx  reaches  an  extreme  value.  Hence, 
a  critical  situation  prevails,  similar  to  the  elastic  case,  if 
0  =  (cp*)l/4  ;  the  corresponding  critical  load  parameter  is 


wr 


On  the  basis  of  the  two  extreme  assumptions  regarding  the  magnitude  of 
F  ,  the  critical  stress  is  assumed  to  be  bounded  by 


✓  E  +  /  Efc 


2  F  EE 

'l  3c  t 

1'+  ,  Et’ 


(125) 


CLASSICAL  MODE 

The  classical  mode  of  instability  is  characterized  by  a  negligible  trans 
verse  strain  and  a  wavelength  that  is  several  times  the  core  thickness. 
The  ratio  of  bending  strains  and  neutral  surface  strains  is  very  small 
for  the  face  layers. 


102 


The  essential  changes  in  the  state  of  strain  are 


U  7r/X  sin  7T* x/X  and 


W/R  sin  7T'x/X 


Hence, 


2  (1  +  p*)  -\/T 


X 

1  +  — 

7rR 


7TX 
sin  — 

X 


Now  as  53  Wg  ,  but  ,  the  term  in  front  of  the  square 

bracket  and  the  second  term  between  the  brackets  have  a  different  sign 
for  the  two  face  layers.  Then  the  sign  of  &eeff  will  also  differ 
for  the  two  faces  if 


7TR 

This  condition  is  met  in  most  practical  cases,  and  then  6cr  will  have 
opposite  signs  in  the  two  face  layers.  The  buckling  criterion  for  the 
classical  mode,  according  to  Equation  (44),  can  be  written  as 


(126) 


and  the  appropriate  expressions  for  the  reduced  bending  and  hoop  stiff¬ 
nesses  of  the  sandwich  shell  are 


-  2 


E  E. 


E  +  E. 


105 


and 


*  , 
E  t 


R 


* 

E  t 


E  +  Et  t 


If  two  bounds  for  the  core  shear  stiffness  are  chosen  as 


G1  ~  <  G1  <  G1 
E 


then  all  stiffness  values  are  stationary  whenever  the  value  of  the  load 
is  stationary  and  bounds  for  the  classical  buckling  stress  are  ob¬ 
tained  as 


c  +  t 


/E  +  E 


1  EE 


E  +  E. 


-i  — p 

ct  E  +  E  E 

- -  - ^  .  —  <  a 

t  2R2  E  +  Et  G1  X 


c  +  t  /E  +  E _  ct  E  +  E  E  E 

- - J  E  Et  -  —  Z—T  — 

E  +  Et  c  2R  E  +  Et  Gx 


(127) 


ANTISYMMETRIC  WRINKLING 


Although,  in  general,  the  critical  stress  for  elastic  antisymmetric 
wrinkling  is  shown  to  be  larger  than  the  elastic  symmetric  wrinkling 
stress,  the  antisymmetric  result  depends  on  both  F  and  G  ,  while 
the  symmetric  result  is  Independent  of  G  ,  Hence,  if  both  quanti¬ 
ties  are  affected  by  plastic  deformations,  the  ratio  of  symmetric  and 
antisymmetric  wrinkling  stresses  will  change.  In  the  antisymmetric¬ 
mode,  face  bending  strains  are  usually  larger  than  face  neutral  sur¬ 
face  strains.  If,  again,  upper  and  lower  bounds  for  the  core  stiffness 
quantities  are  assumed  in  the  form 


®t  .  „  *  *  . 

_  (F ,  G)  '  F  ,  G  <  F  ,  G 

E 
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then  all  stiffness  quantities  assume  stationary  values  at  the  critical 
value  of  the  stress.  Equation  (1*9)  can  be  written  in  the  form 


(128) 


and  the  antisymmetric  wrinkling  stress  is  assumed  to  be  bounded  by 


2(c  +  c)  -y /f 


FEE. 


E  E. 


•v  E  +  E 


-  !*■  2 

t  C  G1 


a  e  +  Et 


(129) 


It  must  be  noted  that  the  prevalence  of  face  bending  strain  over  face 
stretching  is  less  pronounced  in  the  antisymmetric  wrinkling  mode  than 
it  is  in  the  general  instability  mode.  If  the  two  types  of  strain  be¬ 
came  of  comparable  magnitude,  the  change  in  the  stress  intensity  is  a 
more  complicated  function  of  the  coordinates  and  a  more  complete  study 
of  the  Inelastic  behavior  is  mandatory. 
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